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Major  Department:  Physics 

We  present  a  self-consistent  theory  of  superconductors  in  the  presence  of 

Kondo  impurities.  The  impurity  degrees  of  freedom  are  treated  using  the  large- 

N  slave  boson  technique,  leading  to  tractable  equations  describing  the  interplay 

between  the  Kondo  effect  and  superconductivity.    We  show  that  for  a  single 

impurity  in  a  superconductor  with  density  of  states  N(u)  ~  |u;|r,  there  exists 

a  critical  coupling  Jc  below  which  the  Kondo  effect  does  not  occur.  However, 

for  r  <  1  or  N  =  2  any  finite  concentration  of  impurities  drives  Jc  — >  0.  The 

technique  is  tested  on  the  s-wave  case  and  shown  to  give  good  results  compared 

to  other  methods  for  Tj<>Tc.  We  calculate  low  temperature  thermodynamic 
and  transport  properties  for  various  superconducting  states,  including  isotropic 
s-wave  and  representative  anisotropic  model  states  with  line  and  point  nodes  on 
the  Fermi  surface.  The  theory  provides  support  for  phenomenological  models 
of  resonant  impurity  scattering  in  heavy  fermion  systems. 

Motivated  by  some  recent  experiments  on  high-Tc  superconductors,  we 
study  properties  of  superconducting  states  with  "extended  s-wave"  symmetry. 
In  the  presence  of  impurities,  thermodynamic  properties  of  such  states  exhibit 
gap  less  behavior,  reflecting  a  residual  density  of  states.   While  for  a  range  of 

v 


impurity  concentrations,  properties  reflecting  the  density  of  states  alone  will 
be  similar  to  those  of  d-wave  states,  transport  measurements  may  be  shown  to 
qualitatively  distinguish  between  the  two. 

We  also  discuss  the  effect  of  two-channel  Kondo  impurities  on  supercon- 
ductivity. In  the  strong  coupling  regime  such  impurities  are  pairbreakers,  in 
contrast  to  the  ordinary  Kondo  effect.  Measurements  of  Tc-suppression  may 
help  in  identifying  impurities  displaying  this  more  exotic  exchange  coupling  to 
the  conduction  band. 
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CHAPTER  1 
INTRODUCTION 

The  problem  of  a  magnetic  impurity  in  a  superconductor  has  been  exten- 
sively studied,  but  is  not  completely  solved  because  of  the  difficulty  of  treating 
the  dynamical  correlations  of  the  coupled  impurity-conduction  electron  sys- 
tem together  with  pair  correlations.  Generally,  the  behavior  of  the  system 
can  be  characterized  by  the  ratio  of  the  Kondo  energy  scale  in  the  normal 
metal  to  the  superconducting  transition  temperature,  Tj^/Tc.  For  T^/Tc  <C  1, 
conduction  electrons  scatter  from  classical  spins  and  physics  in  this  regime 
can  be  described  by  the  Abrikosov-Gor'kov  theory  [  1].  In  the  opposite  limit, 
Tr'/Tc  3>  1,  the  impurity  spin  is  screened  and  conduction  electrons  undergo 
only  potential  scattering.  In  this  regime  s-wave  superconductors  are  largely 
unaffected  by  the  presence  of  Kondo  impurities  due  to  Anderson's  theorem  [ 
2].  Superconductors  with  an  anisotropic  order  parameter,  e.g.  p-wave,  d-wave 
etc.,  are  strongly  affected,  however,  and  potential  scattering  is  pair-breaking. 
The  effect  of  pair  breaking  is  maximal  in  s-wave  superconductors  in  the  in- 
termediate region,  T%  ~  Tc,  while  in  the  anisotropic  case  it  is  largest  for 
Tr/Tc  ->  co-  The  interest  in  the  latter  class  of  systems  comes  from  studies 
of  the  heavy-fermion  and  high-Tc  superconductors,  exhibiting  many  properties 
that  can  be  explained  by  an  unconventional  order  parameter.  The  word  "un- 
conventional" here  means  that  there  are  broken  symmetries  in  addition  to  the 
U(l)  symmetry  broken  in  classic  superconductor. 
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The  Kondo  effect  is  accompanied  by  the  formation  of  a  narrow  many-body 
resonance  of  width  Tr  near  the  Fermi  level  in  the  impurity  spectral  density 
Nf(io).  It  is  reasonable  to  expect  that  the  opening  of  a  gap  A  in  the  con- 
duction electron  density  of  states  leads  through  hybridization  processes  to  a 
similar  gap  in  NAu),  which  destroys  the  Kondo  effect  if  sufficiently  large. 
Withoff  and  Fradkin  [  3]  pointed  out  that  the  two  problems  of  impurity  spins 
coupled  to  baths  of  conduction  electrons  with  (a)  constant  density  of  states 
and  (b)  with  a  fully  developed  gap  represent  two  extreme  members  of  a  family 
of  problems  given  by  specifying  a  generalized  conduction-electron  density  of 
states  N(u)  =  C|o;|r,  \u\  <  D,  and  0  <  r  <  oo.  Making  use  of  renormaliza- 
tion  group  arguments  as  well  as  explicit  calculations  for  the  large  degeneracy 
SU(iV)  Kondo  (Coqblin-Schrieffer)  model,  they  showed,  for  r  >  0,  the  exis- 
tence of  a  critical  coupling  Jc  below  which  impurities  are  decoupled  from  the 
conduction  band  and  no  Kondo  effect  occurs.  The  potential  physical  examples 
of  this  phenomenon  include  unconventional  superconducting  states  with  line 
and  point  zeroes  in  the  momemtum-dependent  gap  function,  corresponding  to 
densities  of  states  N(oj)  varying  as  u  and  u  ,  respectively.  Such  states  are 
possibly  realized  in  the  heavy  fermion  superconductors  UPt3,  UBei3,  URu2Si2, 
CeCu2Si2,  and  UPd2Al3.  Other  examples  of  systems  which  have  a  gapless  ex- 
citation spectrum  under  certain  conditions  are  a-Sn,  Pbi_xSnxTe  at  a  critical 
composition,  domain  walls  of  PbTe,  PbTe-SnTe  hetero junctions  and  graphite. 

In  the  context  of  heavy-fermion  superconductivity,  theories  of  impurity 
scattering  in  such  states  have  been  given  by  Ueda  and  Rice  [  4]  for  the  case  of 
weak  potential  scattering,  and  by  Hirschfeld,  Vollhardt,  and  Wolfle  [  5],  and 
Schmitt-Rink,  Miyake,  and  Varma  [  6]  for  strong  scattering.  Moment  formation 
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was  not  considered  in  these  theories,  but  pair  breaking  still  occurs  because  of 
the  vanishing  of  the  anomalous  one-electron  impurity-averaged  self-energy.  The 
strength  of  the  scattering  was  parametrized  in  the  latter  works  by  a  phase  shift 
Sq  for  s-wave  potential  scattering  of  electrons  at  the  Fermi  surface.  One  of  the 
principal  results  of  these  treatments  was  that  in  the  resonant  scattering  limit, 
6q  — >  7r/2,  corresponding  to  the  single-impurity  spin-^  Kondo  effect  [  7],  a 
"bound-state"  resonance  was  found  to  form  in  the  superconducting  density  of 
states.,  N(ui),  leading  to  gapless  effects  in  thermodynamic  properties.  These 
are  the  analogs  of  the  bound  states  found  in  discussions  of  Kondo  effects  in 
s-wave  superconductors  [  8]. 

While  this  work  was  begun  in  anticipation  of  applications  to  heavy  fermion 
superconductors,  recent  measurements  of  penetration  depth  [  9],  photoemis- 
sion  [  10,  11],  and  Josephson  tunneling  on  YBCO  and  BSSCO  [  12,  13,  14, 
15]  have  provided  evidence  that  the  copper  oxide  superconductors  may  be  un- 
conventional, possibly  d-wave  as  well.  This  conclusion  remains  controversial, 
however,  and  experimental  tests  to  distinguish  conventional  from  unconven- 
tional pairing  are  of  great  current  interest.  One  would  like,  for  example,  to 
develop  an  understanding  of  the  effect  of  doping  Zn  and  Ni  impurities  in  the 
Cu02  planes.  Simple  models  of  Zn  and  Ni  acting  as  strong  and  weak  potential 
scatterers  in  a  d-wave  superconductor,  respectively,  are  consistent  with  some 
experiments  at  low  doping  levels  [  16,  17,  18],  but  inconsistent  with  other  mea- 
surements [  19].  Since  in  some  cases  Zn  impurities  appear  to  possess  a  magnetic 
moment  at  higher  temperatures  [19],  it  is  of  interest  to  explore  whether  an  s-d 
type  exchange  coupling  of  conduction  electrons  to  an  impurity  embedded  in  an 
unconventional  superconductor  can  describe  the  range  of  behavior  observed. 
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In  this  work  we  focus  initially  on  basic  thermodynamic  and  transport  prop- 
erties of  superconductors  doped  with  Kondo  impurities.  Our  aim  is  to  develop 
a  tractable,  self-consistent  scheme  for  the  calculation  of  all  basic  properties  of 
superconductors  using  methods  known  to  successfully  describe  the  most  diffi- 
cult aspect  of  the  problem,  namely  the  dynamics  of  the  Kondo  impurity.  For 
this  reason  we  have  adopted  the  large- iV  "slave  boson"  approach  of  Barnes 
[  20],  Coleman  [  21],  and  Read  and  Newns  [  22],  as  this  approach  is  well- 
known  to  provide  a  good  description  of  the  spectral  properties  of  the  Kondo 
impurity  for  sufficiently  low  temperatures.  We  emphasize  that  this  description 

is  adequate  for  most  of  our  purposes  {Tc<Tx)  even  in  the  case  of  spin  degen- 
eracy N  =  2,  since  the  impurity  spectral  resonance  is  located  at  the  proper 
position,  i.e.  exactly  at  the  Fermi  level  in  this  approximation. 

Impurity  doping  experiments  in  high  temperature  superconductors  mo- 
tivated us  to  consider  also  the  effect  of  nonmagnetic  impurities  on  strongly 
anisotropic  superconductors  of  the  s-wave  type.  Many  experiments  probe  only 
the  low- lying  quasiparticle  states  of  the  superconductor,  which  in  the  pure  sys- 
tem are  quite  similar  for  both  anisotropic  s-  and  d-wave  states  with  nodes. 
Nonmagnetic  scattering  reduces  anisotropy  of  such  states  without  destroying 
superconductivity,  in  contrast  to  the  effect  the  nonmagnetic  impurities  have 
on  unconventional  (non-s-wave)  superconducting  states.  In  anisotropic  s-wave 
states  having  lines  of  zeroes  of  the  order  parameter  on  the  Fermi  surface,  some- 
times called  "extended-s"  states,  scattering  on  such  impurities  removes  the 
nodes  of  the  order  parameter  at  certain  critical  concentration,  and  induces 
an  energy  gap  at  higher  concentrations.   This  gap  should  be  visible  in  low-T, 
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low  energy  experiments.  If  the  gap  is  too  small  to  be  observed  in  thermody- 
namic experiments,  one  should  try  transport  measurements,  e.g.    microwave 
conductivity  or  thermal  conductivty. 

In  chapter  2  we  address  questions  related  to  the  conditions  for  the  oc- 
curence of  the  Kondo  effect  in  systems  with  the  density  of  states  vanishing  at 
the  Fermi  level  [  23].  In  chapter  3  we  calculate  basic  properties  for  both  con- 
ventional and  unconventional  superconductors  doped  with  Kondo  impurities 
[  24].  Methods  of  distinguishing  strongly  anisotropic  s-wave  superconductors 
from  d-wave  superconductors  having  very  similar  quasiparticle  excitation  spec- 
trum are  discussed  in  chapter  4  [  25,  26].  In  chapter  5  we  study  the  effect 
of  two-channel  Kondo  impurities  on  superconductors  [  27].  Conclusions  are 
presented  in  chapter  6. 


CHAPTER  2 

KONDO  EFFECT  IN  SYSTEMS 

WITH  DENSITY  OF  STATES 

VANISHING  AT  THE  FERMI  ENERGY 


2.1  Mean-Field  Theory  of  the  Single-Impurity  Problem 
An  antiferromagnetic  exchange  interaction  between  a  magnetic  moment  of 
an  impurity  and  a  single  channel  of  conduction  electrons  leads  [  28,  29]  at  low 
temperature  to  the  formation  of  a  spin  singlet  state.  The  temperature  scale 
over  which  the  crossover  to  the  low-T  state  with  compensated  impurity  spin 
occurs  is  called  the  Kondo  temperature  Tr.  The  resistivity  per  impurity  as  a 
function  of  temperature  becomes  logarithmically  divergent  when  the  crossover 
is  reached  and  saturates  at  low  T,  remaining  finite  at  T  =  0.  The  low-T  state 
is  a  Fermi  liquid  with  renormalized  Fermi  liquid  parameters.  [28,29]  It  is  well 
known  that  any  nonzero  conduction  electron  density  of  states  at  the  Fermi 
level  N(ep)  leads  to  the  Fermi  liquid  ground  state  in  the  Kondo  problem. 

What  happens  in  systems  with  N(e)  vanishing  at  e  =  epl  Is  the  ground 
state  of  such  system  containing  a  dilute  ensemble  of  Kondo  impurities  metallic 
or  not?  If  the  low  energy  density  of  states  is  described  by  a  power  law,  N(e)  ~ 
C|e|r,  how  does  the  low-temperature,  low-energy  physics  depend  on  J  and  r? 
Examples  of  systems  in  which  the  density  of  states  vanishes  with  a  power  law 
near  ep  are  unconventional  superconductors  and  semimetals.  To  have  a  basic 
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understanding  of  the  problem,  Withoff  and  Fradkin  [3]  used  poor-man's  scaling 
for  the  spin- ^  Kondo  model 

(2.1) 


+  oJ^Clfk'iS     +4;Cfe/tS+), 


2 

k,k' 

with  the  power-law  density  of  states  over  the  whole  bandwidth.  In  poor-man's 
scaling  [  30]  one  eliminates  high-energy  states  near  band  edges,  reducing  the 
width  of  the  band  to  D  —  SE.  The  new  effective  coupling  constant  Jeyy  is  given 
by  J1  =  J  +  N(  —  D).J'SE/D  +  ...  To  restore  the  original  number  of  states, 
20D1+r/(l  +  r),  the  unit  of  length  is  scaled  by  (D/D')1+r.  The  coupling 
constant  is  changed  by  (D' / D)  +r .  To  restore  the  cutoff  all  quantities  with 
units  of  energy  are  multiplied  by  D/D'.  The  renormalized  coupling  constant 

Jr  is  then 

/ D'V 
Jr  =  f  _  J    J'  ~  J  +  J(JCDr  -  r)SE/D.  (2.2) 

The  fixed  points  are  J  =  0,  J  =  oo,  and  Jc  =  r/(CDr).  The  last  fixed  point 
is  unstable.  For  J  <  Jc  the  effective  J  flows  to  0,  and  at  J  >  Jc  the  coupling 
constant  flows  to  oo.  In  the  ordinary  Kondo  problem  with  N(e)  =  const., 
there  are  only  two  fixed  points:  the  unstable  one  at  J  =  0  and  the  stable  one 
at  J  =  oo.  This  is  then  an  indication  of  qualitatively  different  physics  in  the 
problem  with  r  ^  0. 

Is  the  result  of  the  poor-man's  scaling  confirmed  by  other  approaches,  e.g. 
the  large- iV  approximation?  The  slave  boson  approximation  is  known  to  give 
correct  results  in  the  limit  of  large  orbital  degeneracy  for  the  ordinary  Kondo 
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problem  at  low  temperature  [20,21].    The  simplest  version  of  the  Anderson 
Hamiltonian  to  describe  a  rare-earth  impurity  is 

H  =  Yl  ekclmckm  +  EfYl  f™fm  +  U    S    fmfmfinffmf 

km  m  m>m'  ,~  ~\ 

+  Y,Vk(CLf™  +  h.C.). 
km 

The  electrons  are  labelled  by  their  magnetic  quantum  number  m  =  — j,  ...,j, 
with  j  the  total  angular  momentum,  and  the  impurity  degeneracy  is  N  =  2j  +  l. 
Here  we  ignore  spin-orbit  interaction  and  crystal  fields.  An  electron  on  one  of 
the  degenerate  /  orbitals  at  the  impurity  site  has  energy  Et.  The  energy  U  of 
the  Coulomb  repulsion  between  two  impurity  /-electrons  may  be  large,  e.g.  5- 
10  eV  in  the  lanthanides.  The  impurity  electrons  hybridize  with  the  conduction 
band  with  amplitude  Vk.  Here  we  assume  Vfc  =  V  =  const.,  independent  of 
fe.  When  the  bare  impurity  level  Et  lies  deep  below  the  Fermi  energy,  such 
that  \V2/Ef\  <g.  D,  and  for  large  on-site  Coulomb  repulsion,  U  »  V2/D, 
the  Hamiltonian  (2.3)  can  be  simplified  by  the  Schrieffer- Wolff  transformation 
[  31].  The  resultant  Hamiltonian  is  the  SV(N)  generalization  of  the  Kondo 
model  (2.1),  also  called  the  Coqblin-Schrieffer  model  [  32], 

H  =  J2eAmCkm  +  ^  YY1  C\mfm>frnCkm,.  (2.4) 

km  kk'  mm' 

The  four-fermion  term  can  be  eliminated  by  the  Hubbard-Stratonovich  trans- 
formation. In  the  limit  N  — >  oo,  the  saddle-point  solution  is  equivalent  to  a 
mean  field  theory  with  the  boson  field  playing  the  role  of  an  (unphysical)  order 
parameter,  [22] 


°=jjY,V™ekTn),  (2.5) 

km 
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The  mean  field  Hamiltonian  now  reads 

H  =  Ee*4mc*™  +  E  (°Aj™  +  *-c)  +  ef(nf  -  1).  (2.6) 

km  km 

The  auxiliary  field  and  position  of  the  saddle  point  are  obtained  by  solving  two 

self-consistent  equations 


dH 


=  0,  (2.7) 


and 


8H\ 


'do  J 


=  0.  (2.8) 


For  a  constant  conduction  electron  density  of  states,  N(u)  =  Nq  =  1/2D,  |u;|  < 
D,  equations  (2.7)  and  (2.8)  were  solved  by  Read  and  Newns  [22],  leading  to  a 
Lorentzian  impurity  spectral  density  centered  at  et,  of  width  T  =  ttNqo  .  The 
Kondo  temperature  in  their  approach  is  the  width  of  this  low  energy  resonant 
level,  given  by  TK  =  Jv2  +  e\  =  Dexy>(-1/NN0J).  Equation  (2.7)  can  be 
written  as 

^  =  J^def(e)N(e)Nf(e),  (2.9) 

where  Nj(e)  =  —^lmGj(e+iO+)  is  the  impurity  density  of  states,  and  G7  (u>)  = 
u>  —  €f  —  Sm(cj)  is  /-electron  Green's  function,  while  equation  (2.8)  becomes 


Na 


=   f"  def(e)N(e)- '    1,    v    ,  „2-  (2.10) 
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Figure    2.1     Schematic  flow  diagram  for  single-channel  finite- r  Kondo  prob- 
lem obtained  from  the  Numerical  Renormalization  Group  (based  on  Ref.  [26] ) . 


The  impurity  self-energy  Sm 


is 


.W^I^l-taAvW. 


(2.11) 


The  critical  coupling  obtained  from  the  solution  of  equations  (2.9)  and  (2.10) 
corresponds  to  a  =  e  /  =  0  at  T  =  0.  From  equation  (2.10)  the  critical  coupling 
is  Jc  ~  r/CDr .  Both  a  and  ey  vanish  as  power-law  functions  of  \J  —  Jc\  near 
Jc,  a  ~  (J  —  Jc)1'r,ey  ~  (J  —  Jc)l'r-  This  result  comes  from  the  work  of 
Withoff  and  Fradkin  and  was  later  confirmed  by  us  in  a  numerical  solution  of 
equations  (2.9)  and  (2.10).  The  energy  scale  for  the  formation  of  the  singlet  is 
given  near  the  transition  by  Tr  =  (r(J  —  Jc)/CJJc)'r). 

Recent  Numerical  Renormalization  Group  calculations  confirm  the  mean- 
field  result  for  r  <§C  1  but  differ  substantially  at  larger  r  [33].  The  qualitative 
form  of  the  flow  diagram  obtained  within  NRG  is  shown  in  Figure  2.1  .  There 
is  a  critical  exponent,  rc  =  0.5,  such  that  for  r  >  rc  the  Kondo  effect  does  not 


occur. 
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In  a  more  realistic  treatment,  the  density  of  states  may  be  assumed  to  vary 
as  a  power-law  over  a  small  energy  scale  A  <C  D,  N(u>)  =  C\u\r,  and  to  be 
constant  otherwise,  N(co)  =  C,  for  A  <  |cu|  <  D.  We  find  that  the  critical 
coupling  obtained  in  this  case  is  independent  of  r,  at  least  for  r  >  1,  Jc  — 
2D/\n(2D/A).  It  was  shown  for  the  problem  with  a  small  gap,  N(lj)  =  0  for 
\u)\  <  A,  that  Quantum  Monte  Carlo  and  Numerical  Renormalization  Group 
approach  yield  similar  results  [  34,  35,  36,  37,  38].  The  transition  associated 
with  the  critical  ratio  T^/A  might  be  observable  provided  we  can  tune  T/^/A. 
However,  it  is  not  a  phase  transition  in  the  usual  sense  -  it  can  occur  only 
at  T  =  0.  At  any  finite  T  the  system  will  not  exhibit  singularities  of  the 
thermodynamic  functions  but  it  will  have  instead  a  smooth  crossover  from 
the  impurity  spin-singlet-like  characteristics  to  a  free  moment  behavior  with 
increasing  A.  One  of  the  quantities  of  interest  to  experimental  studies  of  such 
a  transition  is  the  spin  susceptibility.  In  the  ordinary,  single  channel  Kondo 
problem,  the  static  susceptibility  reaches  the  limit  \  — >  const. /Tj^,  or  T\  — >•  0, 
as  T  — >  0,  indicating  the  screening  of  the  impurity  spin  by  conduction  electrons. 
In  a  system  with  a  gap  in  the  conduction  electron  spectrum,  the  formation  of 
a  singlet  state  is  less  advantageous  energetically  and  \(T  =  0)  becomes  finite 
as  the  gap  size  increases  or  electron  density  decreases.  Measurements  of  \  at 
finite  T  should  indicate  a  crossover  to  a  Curie-like  suceptibility  as  the  singlet 
ground  state  becomes  unstable. 

2.2  NCA  Approximation 
A  more  sophisticated  version  of  the  large-TV  approach,  the  noncrossing 
approximation  (for  review  of  the  method  see  Refs.      [  39]  and    [  40])  leads  to 
the  same  conclusions  and  allows  for  an  extension  of  the  theory  into  the  high 
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Nf(w) 


Lf  fj. 

Figure  2.2  Schematic  impurity  density  of  states  for  the  Anderson  model;  fx 
is  the  Fermi  level  and  Ej  is  the  position  of  the  bare  impurity  level.  The  Kondo 
temperature  is  the  width  of  the  low  energy  many-body  resonance. 


(2.12) 


temperature  regime,  T  S>  T^.  In  NCA  one  sums  all  diagrams  in  the  projected 

perturbation  expansion  with  no  crossing  conduction  electron  lines.  Here  we  do 

the  calculations  for  the  Anderson  model, 

H  =  Yl  ekclmckm  +  EfYl  fmfm  +  VY1  [cL/™ftt  +  h'c- 
k,m  m  k,m 

where  A  is  a  Lagrange  multiplier  enforcing  the  constraint  nj+n^  =  Ylm  fmfm  + 
b+b  =  1,  preventing  double  occupancy  of  the  impurity  site. 

The  impurity  density  of  states  for  the  Anderson  model  in  the  Kondo  regime, 
— Ef  ^>  V  /D  is  shown  schematically  in  Figure  2.2  .  The  narrow  feature 
at  the  Fermi  level  is  the  Abrikosov-Suhl  many-body  resonance  of  width  Tr, 
which  controls  all  low-T  thermodynamic  and  transport  properties.  The  NCA 
is  known  to  reproduce  all  qualitative  features  of  the  Anderson  model  spectral 
function  on  the  scale  shown,  and  fails  only  below  an  energy  T^ca  ~  T^/D  [ 
41,40]. 

The  self-energies  for  the  boson  and  fermion  propagators  are,  respectively 

(see  Figure  2.3  ), 

/oo 
def{e)N(e)Gm(oj  +  e  +  iO+),  (2.13) 

-00 
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Figure    2.3     Self-energies  for  the  empty  impurity  /-level  (slave  boson),  Sq, 
and  for  the  occupied  /-level  (fermion),  Hm. 


and 

/oo 
de(l  -  f(e))N(e)G0(u  -  e  +  i0+).  (2.14) 

-oo 

At  T  =  0  and  for  integer  r  equations  (2.13)  and  (2.14)  can  be  tranformed  into 
a  system  of  differential  equations  and  can  be  solved  by  a  generalization  of  the 
method  used  by  Inagaki  [  42]  and  later  by  Kuramoto  and  Kojima  [  43].  The 
solution  of  NCA  equations  at  finite  temperature  for  r  =  0  is  discussed  in  detail 
e.g.  in  a  review  by  Bickers  [40]. 

Here  we  only  show  the  results  for  the  spin  susceptibility.  The  temperature 
dependence  of  the  static  spin  susceptibility  exhibits  a  crossover  from  the  high-T 
free-moment  behavior  to  the  low-T  strong  coupling  regime. 

The  levelling  off  of  the  r  =  0  curve  in  Figure  2.4  indicates  a  crossover  to 
the  low  temperature  Fermi  liquid  state,  as  confirmed  by  the  Bethe  Ansatz  [  44] 
and  NRG  solutions  of  the  problem  [29].  Figure  2.4  captures  the  essence  of  the 
problem  correctly:  increasing  r  leads  to  lower  Tg . 
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Figure  2.4  The  spin  susceptibility  for  r  =  0  (circles),  0.1  (triangles),  0.2 
(rhombuses);  N  =  2,  T/D  =  0.25,  Ef/D  =  -0.67.  The  Kondo  tempera- 
ture is  TK/D  =  (r/7r£>)1/^exp(7rE//iVr),  and  for  r  =  0,  T/D  =  0.25,  it  i 
approximately  equal  to  Tr/D  ~  4.2  x  10-3. 


is 


In  Figures  2.5  ,  2.6  ,  2.7  ,  and  2.8  we  show  log\(T,u;  =  0)  for  r  =  0.05,0.1, 
0.15,  0.2,  and  several  values  of  hybridization.  We  clearly  see  a  decrease  of  T^ 
with  decreasing  V  in  each  set  of  data. 

Figure  2.9  shows  log(T\)  for  two  lowest  values  of  T  from  each  of  the  Figures 
2.5-2.8.  In  the  case  of  r  =  0.2  the  transition  from  the  screened  to  the  finite 
moment  regime  is  rather  convincingly  demonstrated.  For  r  =  0.05,0.1,  and 
0.15,  the  values  of  T  seem  not  to  be  far  from  the  critical  ones,  although  the 
transition  is  not  reached.  How  does  Fc  obtained  in  NCA  compare  to  the 
critical  coupling  obtained  in  the  poor  man's  scaling  and  in  the  mean  field 
approach?  Remembering  that  in  the  CS  limit  J  =  V2/Et  =  V/nEt,  we  find 
that  iVVc,NCA(r  =  0.2)  ~  0.16,  where  we  assumed  Tc  ~  0.32  for  r  =  0.2. 
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Figure  2.5  The  spin  suceptibility  for  r  =  0  and  r  =  0.05;  TV  =  2,  Ef/D  = 
-0.67.  For  r  =  0  and  T  =  0.14,0.15,0.175,  and  0.2,  the  Kondo  temperature  is 
TK/D  ~  1.15  x  10~4, 1.96  x  10-4,5.77  x  10~4,  and  1.31  x  10"3,  respectively. 


The  critical  coupling  N0JC  for  r  =  0.2,  in  both  the  mean  field  and  the  NRG 
approach  is  approximately  0.2. 
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Figure  2.6  The  spin  suceptibility  for  r  =  0  and  r  =  0.1;  N  -  2,  £y/.D  = 
-0.67.  For  r  =  0  and  r/D  =  0.2,0.225,0.25,  and  0.275,  the  Kondo  tem- 
perature is  TK/D  ~  1.31  x  10~3,2.5  x  10~3,4.2  x  10~3,  and  6.44  x  10-3, 
respectively. 
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Figure    2.7     The  spin  suceptibility  for  r  =  0  and  r  =  0.15;  N  —  2,  Et/D  = 

-0.67.  For  r  =  0  and  r/D  =  0.3,  0.35,  0.4,  0.5,  and  0.6,  the  Kondo  temperature 

is  TK/D  =  9.3  x  10-3,1.65  x  10~2,2.6  x  10~2,4.9  x  10~2,  and  7.6  x  10~2, 
respectively. 
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Figure  2.8  The  spin  suceptibility  for  r  =  0  and  r  =  0.2;  iV  =  2,  £"//-D  = 
-0.67.  For  r  =  0  and  T/D  =  0.3,  0.35,  0.4,  and  0.5,  the  Kondo  temperature  is 
TK/D  =  9.3  x  10"3, 1.65  x  10_2,2.6  x  10"2,  and  4.9  x  10-2,  respectively. 
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Figure   2.9     Logarithm  of  T\  for  two  lowest  values  of  T  from  each  set  of  data 
shown  in  Figures  2.5-2.8. 


2.3  Kondo  Impurity  in  a  Superconductor 
One  of  the  obvious  applications  of  the  theory  discussed  in  the  previous 
section  is  the  problem  of  a  Kondo  impurity  in  a.  superconductor.  A  supercon- 
ductor differs  from  a  normal  metal  not  only  through  the  gap  in  the  density  of 
states,  but  also  through  the  existence  of  off-diagonal  long  range  order  described 
by  the  order  parameter  A(k).  To  analyze  the  problem  we  start  from  the  SU(N) 
Coqblin-Schrieffer  Hamiltonian,  already  introduced  in  equation  (2.4),  and  add 
a  simple  BCS-like  pairing  of  electrons  on  opposite  sides  of  the  Fermi  sphere, 
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H  =  Yl  eAmCkm  +xY,I2  ClmfL'f^km'  +  Y,  [A(kHmC-k-m  +  hx- 
k,m  k,k'  m,m'  k,m 

(2-15) 

We  now  generalize  the  procedure  of  Read  and  Newns  [22]  to  include  super- 
conducting correlations  in  the  functional  integral  representation  of  (2.15).  The 
saddle-point  approximation  to  this  theory  is  equivalent  in  the  N  — >•  oo  limit  to  a 
mean-field  theory  of  (2.15)  with  mean-field  amplitude  a  =  (J/N)  Ylkm(cmfm) 
and  Lagrange  multiplier  et  implementing  the  average  constraint.  It  leads  to 
the  two  equations 


f°°  1 

_  =  -Imy_     duf(u)-Tr[(T0  +  T3)Gf(u  +  iO+)],  (2.16) 

and 


/°°  1       r  i 

dujf{uj)-Tt   (T0  +  TS)G°(k,uj  +  iO+)Gf(uj  +  iO+)    .       (2.17) 
-oo  z        L  J 


where  G°  denotes  the  conduction  electron  Green's  function  in  the  pure  super- 
conductor and  Gy  is  the  full  impurity  Green's  function,  given  by 

GJ1(uj)  =  uT0-efT3--Ef(uj).  (2.18) 

Both  are  matrices  in  particle-hole  space  spanned  by  the  Pauli  matrices  rz-, 

G=(#     £)      and     G/=(^t     ^  (2.19) 

and  Sj(u;)  =  a  z2kG(k,u)  is  the  impurity  self-energy.  In  the  supercon- 
ducting state,  we  must  solve  the  full  saddle-point  equations  (2.13)  and  (2.14) 
together  with  Dyson's  equations  for  G y  and  G  as  well  as  the  gap  equation 
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/oo 
duf(u)  Y,  Vk'Tr-  [(n  -  ir2)G(k',  u)]  .  (2.20) 

k'  Z 

Let  us  focus  first  on  the  case  of  a  single  impurity.    We  also  consider  the 

case  of  small  N,  despite  the  fact  that  Eqs  (2)  and  (3)  are  strictly  valid  only  in 
the  large- iV  limit.  Here  we  adopt  the  point  of  view  that,  since  the  saddle  point 
for  N  =  2  is  known  to  reproduce  the  correct  analytic  low-temperature  normal- 
state  behavior  [7]  of  the  /  resonance,  including  its  position  at  €f  =  0,  the 
N  =  2  theory  will  provide  a  good  starting  point  for  a  description  of  Tr  3>  T 
regime.  For  the  superconducting  order  parameter  we  take  for  simplicity  the 
usual  isotropic  s-wave  state  A(fc)  =  A,  model  p-wave  states  with  lines  ["polar", 
A(k)  =  Ao&z]  and  points  of  nodes  ["axial",  A(fc)  =  &o(kx+iky)]  on  the  Fermi 
surface,  with  densities  of  states  varying  at  low  energies  u>  <C  Aq  as  0,  uj  and  u)  , 
respectively.  The  critical  coupling  Jc  in  this  case  [3]  is  now  defined  to  be  that 
J  for  which  a  =  ej  =  0  is  the  only  solution  of  Eqs.  (2)  and  (3),  with  G(k,ui) 
replaced  by  G  (k,u).  We  note  that  to  show  that  a  and  et  always  scale  to  zero 
together  at  the  transition  for  N  <  oo  requires  a  careful  analysis  of  impurity 
bound  states  in  Af(u),  which  occur  in  the  gap  and  outside  the  band  edges.  It 
follows  from  this  analysis  and  from  Eq.  (3)  that  Jc  is  independent  of  N.  For  the 
(unphysical)  case  Aq  =  D  we  recover  the  WF  r  =  1  result,  Jc/D  =  1,  for  the 
polar  state,  while  for  the  axial  state  we  find  Jc/D  =  1.44.  This  differs  slightly 
from  the  WF  r  =  2  result  Jc/D  =  1.33,  as  the  axial  density  of  states  deviates 
from  pure  u  behavior  at  larger  energies.  In  the  physical  limit  Aq  <C  D,  we 
obtain  Jc  ~  2D/\ii(2D/Aq)  for  axial,  polar  and  isotropic  s-wave  states. 
In  Figure  2.10  we  plot  Jc  vs  Aq/D  for  all  three  states. 
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Figure     2.10     Critical  coupling  Jc/D  for  one  impurity  as  a  function  of  the 
order  parameter  amplitude. 


To  study  the  case  of  finite  impurity  concentrations,  we  calculate  self-consistent 
Green's  functions  averaged  over  impurity  positions  in  the  usual  way,  leading 
to  G7  (to)  =  lotq  —  e/T3,  and  G      (u;)  =  ujtq  —  e^r^  —  A(k)r\,  where 


Q  =  u  +  r  lib  I  (A2(k)  -  lD2)    2\  _ 


and 


u 


-.2 


u)  +  aTuj/{— u    +€f). 


(2.21) 


(2.22) 


Here  a  =  nT^N /2n,  and  n  =  7i/TcqNq  is  the  scaled  impurity  concentration. 
In  general,  G7  and  G  will  also  contain  additional  off-diagonal  renormal- 
izations,  which  vanish  in  the  p-wave  case  considered  here. 

An  interesting  consequence  of  the  self-consistent  treatment  of  impurity 
scattering  is  that  the  conclusions  of  WF  regarding  the  existence  of  a  critical 
exchange  coupling  Jc  based  on  a  single-impurity  analysis  are  modified.  It  is 
clear  from  physical  considerations  or  from  Eqs.  (2)  and  (3)  that  the  Kondo 
effect  occurs  for  all  J  >  0  whenever  the  density  of  states  at  the  Fermi  level 
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N(0)  is  finite.  Since  in  the  polar  state  any  finite  impurity  concentration  n  may 
be  shown  to  lead  self-consistently  to  N(0)  >  0,  as  also  found  in  Refs.  [4,5,6], 
the  transition  discussed  by  WF  does  not  take  place.  This  may  easily  be  seen 
by  solving  the  equations  for  <I>  and  u  at  u  =  0,  with  ImiD(O)  >  0.  A  closer 
analysis  shows  that  Jc  =  0  for  all  superconducting  states  with  density  of  states 
N(u>)  ~  |u;|r  for  u  -C  Aq,  r  <  1.  In  the  axial  state  (r  =  2),  and  indeed  for 
any  state  with  1  <  r  <  oo,  the  cases  N  =  2  and  N  >  2  are  qualitatively 
different.  If  N  >  2,  a  critical  concentration  is  required  to  create  a  gapless  state 
N(0)  >  0,  and  thus  drive  Jc  ->  0.  However,  when  the  bound  state  is  located 
exactly  at  the  Fermi  surface  (N  =  2,ey  =  0),  we  find  again  a  finite  density  of 
states  N(0)  for  any  finite  concentration.  These  results  are  also  in  accord  with 
earlier  studies  [5,6],  where  the  phenomenological  phase  shift  8q  is  crudely  given 
here  by  cot-1ey/T. 

Our  results  suggest  that  the  transition  discussed  by  WF  might  be  ob- 
servable in  ordinary  superconductors  doped  with  Kondo  impurities  with  spin 
degeneracy  N  >  2,  e.g.,  Ce.  For  N  =  2  we  have  shown  that  effective  Kondo 
temperature  in  the  superconducting  state  is  reduced  but  never  vanishes.  Never- 
theless, in  relatively  clean  systems  deviations  in  thermodynamic  properties  per 
impurity  from  those  of  the  pure  superconductor  may  be  qualitatively  similar 
to  what  one  might  expect  from  a  WF-type  analysis  if  the  effective  TK  is  driven 
to  zero.  The  theory  presented  here  provides  an  easily  tractable  framework  to 
calculate  such  properties,  as  well  as  transport  coefficients,  in  the  superconduct- 
ing state.  Furthermore,  it  improves  upon  phenomenological  theories  [5,6]  by 
including  a  Kondo  impurity  description  of  the  energy  dependence  of  scattering 
phase  shifts. 


CHAPTER  3 

PHYSICAL  PROPERTIES  OF  MAGNETIC 

IMPURITIES  IN  SUPERCONDUCTORS 


3.1  Introduction 
In  this  chapter  we  focus  on  basic  thermodynamic  and  transport  properties 
of  superconductors  doped  with  Kondo  impurities.  We  would  like  to  examine 
these  properties  both  near  Tc  and  at  T  <C  Tc,  when  the  superconducting  gap  is 
fully  developed.  The  presence  of  bound  states  within  the  superconducting  gap 
strongly  modifies  the  low-T  behavior.  This  is  true  for  both  conventional  and 
unconventional  superconductors.  After  testing  our  approach  and  providing 
some  new  results  for  the  s-wave  case  we  discuss  results  for  unconventional 
superconductors.   We  emphasize  that  the  large-iV  description  is  adequate  for 

most  of  our  purposes  {Tc<Tk)  even  in  the  case  of  spin  degeneracy  N  =  2,  since 
the  impurity  spectral  resonance  is  located  at  the  proper  position,  i.e.  exactly 
at  the  Fermi  level  in  this  approximation. 

As  in  chapter  2,  we  use  the  large- A^  slave  boson  technique  for  the  SU(N) 
Anderson  model  describing  an  N-fold  degenerate  band  of  conduction  electrons, 
Cf.m,  m  =  1,  ...JV  with  energy  e^  hybridizing  through  matrix  element  V  with 
a  localized  impurity  state  fm.  In  general  this  Hamiltonian  contains  a  term 
with  the  Coulomb  repulsion  U  between  two  electrons  present  at  the  impurity 
simultaneously.  In  many  compounds  U  is  large  and  for  the  purpose  of  studying 
the  low  temperature  physics  we  will  assume  U  =  oo.  The  on-site  repulsion  term 

24 


25 
is  then  absent,  but  a  constraint  is  added  to  ensure  that  the  system  remains 
in  the  physical  part  of  the  Hilbert  space.  The  conduction  band  is  assumed  to 
have  a  constant  density  of  states  in  its  normal  state,  N(u)  =  1/2D  =  Nq.  We 
also  include  a  BCS-like  pairing  term  of  electrons  on  opposite  sides  of  the  Fermi 
sphere, 

H  =  H  ek^kmCkm  +  EfYl  f^mfm  +  V  ^[c1 ' kmfmb  +  h.c] 

k.m  m  k,7n 

(3  1) 

+  $}A(*Oct*mCt-*-m  +  h.c]  +  A(]T/tm/m  +  tfb), 

k,m  rn 

In  the  limit  Et  — >  —  oo,  NqV  /Et  =const,  Equation  (3.1)  reduces  to  the 
Coqblin-Schrieffer  Hamiltonian  with  pairing  studied  in  section  2.3.  Here  we 
have  chosen  the  more  general  form  (3.1)  to  study  deviations  from  single  oc- 
cupancy, nt  ^  1,  although  we  do  not  attempt  to  explore  the  fully  developed 
mixed  valent  regime. 

The  mean  field  approximation  to  this  model,  with  mean-field  amplitude 
(6),  leads  to  the  two  equations 


/•oo 

-  =  -lml      du;f(uj)-Tr(T0  +  T3)Gf(u;  +  iO+),  (3.2) 

and 


E  t  —  €  t  Z-00  1         r  i 

fv2  J  =  Imj       duf  (u)-Tr  [(rQ  +  T3)(G°(k, u>  +  iQ+)G}{u  +  iO+))J  , 

(3.3) 
which  determine  (b)  and  et,  the  latter  being  the  position  of  the  resonant  state. 
Equations  (3.2)  and  (3.3)  should  be  solved  self-consistently  together  with  the 
gap  equation  (2.20).    The  full  conduction  electron  Green's  function  G  and 
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Figure  3.1     Dyson  equations  for  the  conduction  electron  and  impurity  Green's 
functions.  Double  lines  are  full  Green's  functions;  single  lines  are  bare  Green's 
functions  (V  =  0). 


the  impurity  Green's  function  Gy  are  calculated  from  the  diagrams  shown  in 
Figure  3.1  ,  yielding 

CM"1  =  G°(ufl  -  E(w)  =  ut0  -  ekr3  -  A(fc)n,  (3.4) 

Gf(u)-1  =  Gfaf1  -  Vf(u)  =  Cotq  -  efr3  -  An-  (3-5) 

The  Green's  functions  are  now  averaged  over  impurity  positions  in  the  usual 

way.    The  renormalized  frequencies  are  calculated  self-consistently  from  the 
Dyson  equations, 


u>  =  u  +  au/{-u2  +  e2f  +  A2),  (3.6) 


and 
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u>  =  u>  +  r/u/(A2{k)-u2y     \.  (3.7) 

Here  a  =  NnTTcQ/2n  ,  n  =  n/N()TcQ  is  the  scaled  impurity  concentration,  and 
(...)r  is  a  Fermi  surface  average.  The  off-diagonal  renormalizations  are 


A{k)  =  A(k)  +  aA/(-u2  +  e2f  +  A2),  (3.8) 


and 


A  =  r/ A(k)/  (A2{k)  -  w2)  1/2\   .  (3.9) 

In  superconductors  with  order  parameters  where  the  Fermi  surface  average 
in  equation  (3.9)  vanishes,  off-diagonal  corrections  vanish  and  A(k)  =  A(k). 
This  class  includes  but  is  not  limited  to  odd-parity  superconducting  states.  The 
energy  scale  T  is  the  renormalized  resonance  width,  T  =  (b)irNoV2.  The  low- 
temperature  Kondo  scale  in  the  large- N  slave  boson  theory  is  given  by  Tj^  = 


T2  +  e/2.  Although  the  width  of  the  actual  spectral  feature  corresponding 
to  the  Abrikosov-Suhl  resonance  is  modified  below  Tc,  in  what  follows  we  will 
normalize  all  quantities  with  respect  to  this  T^-,  evaluated  from  equations  (3.2) 
and  (3.3)  with  A  =  0  at  T  =  0.  In  the  regime  of  principal  interest,  Tjf  >  Tc, 
corrections  to  this  definition  are  small  in  any  case. 

The  early  history  of  the  problem  of  a  Kondo  impurity  in  an  s-wave  su- 
perconductor has  been  reviewed  by  Muller-Hartmann  [  45].  Abrikosov  and 
Gor'kov  first  discussed  the  pair-breaking  effects  of  magnetic  impurities  weakly 
coupled  via  exchange  interactions  to  conduction  electrons  [1].  Shiba  [  46]  ex- 
tended this  approach  to  treat  strong  scattering  by  classical  spins,  using  the 
f-matrix  approximation,  and  showed  the  existence  of  bound  states  in  the  gap. 
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At  finite  concentration  of  impurities  the  bound  states  were  found  to  form  an 
impurity  band  whose  width  and  center  scaled  with  impurity  concentration  and 
exchange  strength.  These  early  works  neglected  the  dynamical  screening  of  the 
localized  spin  by  the  conduction  electron  gas.  These  effects  were  incorporated 
by  Miiller-Hartmann  and  Zittartz  [  47],  adopting  an  equation  of  motion  decou- 
pling scheme  previously  used  by  Nagaoka  [  48]  to  calculate  the  dynamical  spin 
correlations  in  the  normal  state.  This  approach  correctly  reproduced  results 
in  the  Abrikosov-Gor'kov  limit,  Tft/Tc  — >  0,  and  made  the  remarkable  predic- 
tion of  a  reentrant  superconducting  phase  if  Tj(  <C  Tc,  subsequently  observed 
in  experiments  on  Lai_xCexAl2  [  49,  50,  51].  The  failure  of  the  decoupling 
scheme  used  to  capture  the  correct  crossover  to  Fermi  liquid  behavior  in  the 
normal  state  as  T  — >  0  invalidated  the  Muller-Hartmann-Zittartz  approach  in 
the  low  temperature  regime  Tr  >  Tc,  however. 

The  physics  of  the  Fermi  liquid  regime,  T^/Tc  — >•  oo,  was  studied  by 
Matsuura,  Ichinose  and  Nagaoka  [  52]  and  by  Sakurai  [  53]  by  extending  the 
Yamada-Yosida  theory  [  54]  to  the  superconducting  state.  They  obtained  an 
exponential  Tc-suppresssion  with  increasing  impurity  concentration  n,  Tc  ~ 
Tcoexp(—  pn/X),  where  A  is  the  BCS  dimensionless  coupling  constant,  and  p  is 
a  constant  of  order  unity.  This  is  commonly  referred  to  as  "pair- weakening"  as 
opposed  to  pair-breaking,  since  the  effective  superconducting  coupling  constant 
is  reduced  due  to  correlations  on  the  impurity  site.  The  exponential  form 
breaks  down  for  concentrations  sufficiently  close  to  a  critical  nc,  for  which 
Tc  =  0.  In  this  regime  the  reentrant  behavior  found  by  Miiller-Hartmann  and 
Zittartz  does  not  occur.  A  further  characteristic  signature  of  the  Fermi  liquid 
regime  is  the  reduced  specific  heat  jump,  C*  =  (AC/ACo)/(T,c/rco))iji  =7- 
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which  is  always  less  than  one  [53,  55],  in  contrast  to  the  high  temperature 
regime. 

Not  surprisingly,  qualitatively  similar  results  were  obtained  by  other  early 
workers  for  Kondo  and  Anderson  impurities  using  a  variety  of  other  approaches 
[45].  More  recent  treatments  include  the  use  of  a  self-consistent  large- iV  [  56], 
Monte  Carlo  [  57]  and  NRG  methods  [  58].  Schlottman  [  59]  has  treated 
the  mixed-valence  regime  using  Brillouin-Wigner  perturbation  theory.  Out  of 
these  efforts  has  evolved  a  qualitatively  consistent  picture  of  the  effect  of  Kondo 
impurities  on  the  superconducting  transition  [  60],  but  little  understanding  of 
the  low-temperature  properties  of  Kondo-doped  superconductors  because  of 
the  difficulty  of  the  calculations  involved.  In  the  next  section  we  show  that 
the  current  theory  reproduces  the  known  effects  of  Kondo  impurities  on  the 
critical  temperature,  specific  heat  jump  and  bound  states  spectrum  of  an  s- 
wave  superconductor. 

3.2  Results  for  Conventional  Cuperconductors 

3.2.1  Suppression  of  the  Critical  Temperature 

The  simplest  and  most  direct  effect  of  impurity  scattering  on  a  supercon- 
ductor is  the  suppression  of  the  critical  temperature.  Scattering  from  impuri- 
ties with  internal  quantum-mechanical  degrees  of  freedom  leads  to  deviations 
from  the  classic  Abrikosov-Gor'kov  prediction  for  the  dependence  of  Tc  on  im- 
purity concentration  [1].  These  effects  depend  sensitively  on  the  low-energy 
behavior  of  the  self-energy  E(o>),  which  enters  the  linearized  gap  equation, 
obtained  from  equation  (2.20)  near  Tc, 


30 


1.0 


0.5 


0.0 


0.0        20.0       40.0       60.0       80.0      100.0 
n/N0Tc0 


Figure  3.2  The  critical  temperature  for  an  s-wave  superconductor  as  a  func- 
tion of  impurity  concentration.  The  inset  shows  the  slope  of  this  dependence 
evaluated  at  Tc  =  Tcq. 


ln(rc/rc0)  =  2nTc  ]T  * 

^0^(1  +  Q/5W) 


where  B(un)  =  (ujn  +  T)2  +  ef2  . 

The  slope  of  the  Tc-suppression  evaluated  at  ft  =  0  is  therefore 


J_dTc 
Trc\  dn  )  =_ 


N 


E 


rrc0 


n=0  27r  ^  («  +  l/2)S(w„) 


(3.10) 


(3-11) 


In  an  insert  to  Figure  3.2  ,  we  have  plotted  a  numerical  evaluation  of 
equation  (3.11)  for  an  s-wave  superconductor.  Note  the  curve  is  drawn  with  a 
broken  line  for  small  Tj^/Tcq  to  reflect  the  fact  that  the  slave  boson  mean  field 
theory  is  expected  to  break  down  there.  The  maximum  of  the  Tc-suppression 
is  found  to  occur  at  Tr  ~  3TC,  similar  to  the  NCA  result  [56].  The  early 
high-temperature  theory  of  Muller-Hartmann  and  Zittartz  [47]  locates  this 
maximum  at  Tj^jTc  ~  12,  whereas  in  the  more  recent  Monte  Carlo  calculation 
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[57]  for  the  symmetric  Anderson  model  with  finite  U  the  maximum  slope  of  the 
Tc-suppression  is  at  Tj{  ~  Tc.  Unfortunately  a  direct  quantitative  comparison 
with  the  latter  work  is  not  possible,  as  the  simulation  is  not  performed  in  the 
fully  developed  Kondo  regime. 

The  present  theory  predicts  an  exponential  decrease  of  Tc  at  small  concen- 
trations, in  agreement  with  other  theories  of  the  Fermi  liquid  regime  [53,55], 


Tc  ^  Tc0exp 


*  r*(TRr/27rrc0) 


~Tc(l-^HTK/Tco)y        (3-12) 


Tk2 

where  \I>  is  the  digamma  function.  In  Refs.  [  61]  and  [53]  the  initial  suppres- 
sion of  Tc  is  proportional  to  In  (Tj{/Tcq).  A  full  evaluation  of  equation  (3.11) 
for  arbitrary  concentrations  and  various  values  of  the  ratio  Tj^/Tc  is  shown  in 
Figure  3.2  .  It  is  interesting  to  note  that  the  theory  reproduces  the  reentrant 
behavior  characteristic  of  the  high  temperature  regime  [45]  although  we  do  not 
expect  the  theory  to  be  accurate  in  this  case  (dashed  line). 

As  is  evident  in  the  Figure  the  current  theory  predicts  no  critical  concen- 
tration nc  for  which  Tc  =  0.  This  is  a  subtle  point  discussed  by  Sakurai  [53], 
who  suggests  that  a  failure  to  include  the  dynamics  of  magnetic  scattering  by 
states  close  to  the  Fermi  surface  can  lead  to  such  an  effect.  Such  processes  are 
included  in  the  finite-^/  perturbation  theory  through  Coulomb  vertex  correc- 
tions to  the  impurity  averaged  pair  correlation  function.  In  our  U  =  oo  theory, 
such  vertex  corrections  arise  first  in  leading  order  1/iV  corrections  due  to  the 
exchange  of  slave  bosons,  whose  dynamics  are  neglected  in  this  work.  We  ex- 
pect that  effects  arising  from  the  absence  of  these  fluctuations  in  the  theory 

will  be  quantitatively  small  for  Tc<Tj^,  except  for  impurity  concentrations  so 
large  such  that  Tc  <C  Tcq. 
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3.2.2  Position  of  Bound  States 

For  conventional  superconductors,  A(fc)  =  A,  equations  (3.7)  and  (3.9)  are 
simply 


u  =  u  +  T — ,  (3.13) 

(A2-^)1/2 

and 


It  is  obvious  from  equation  (3.14)  that  a  finite  gap  in  the  conduction  electron 
spectrum  induces  a  gap  in  the  impurity  spectrum,  Nf(ui)  =  —  ^ImGf(cj  +  iO+). 
When  the  Abrikosov-Suhl  resonance,  which  develops  at  temperatures  below 
T/<-,  falls  in  the  superconducting  gap,  bound  states  appear  in  the  conduction 
electron  spectrum.  These  peaks  are  placed  symmetrically  relative  to  the  cen- 
ter of  the  gap  and  their  spectral  weight  depends  on  Tj^/Tq  and  on  impurity 
concentration. 

The  density  of  states  for  TV  =  2  and  N  =  4,  exhibiting  pronounced  peaks 
at  the  bound  state  positions,  is  shown  in  Figure  3.3  .  Figure  3.4  shows  N(u>) 
in  the  Coqblin-Schrieffer  limit  for  N  =  2  and  several  values  of  T.  In  the 
high-temperature  regime  the  bound  states  emerge  from  the  edges  of  the  gap 
and  move  towards  the  center  of  the  gap  as  Tj^/Tc  increases  [  62].  In  the  low- 
temperature  limit,  the  bound  states  disappear  into  the  gap  edges  again  [61]. 
Recently,  Shiba  et  al.  [58]  studied  the  position  of  bound  states  in  the  gap  in 
an  s-wave  superconductor  using  the  numerical  renormalization  group  (NRG). 
Their  results  cover  both  the  high-  and  low-temperature  limits  reliably.    The 
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Figure  3.3  Conduction  electron  and  impurity  spectral  functions  in  the 
Coqblin-Schrieffer  limit  in  an  s-wave  superconductor  for  N  =  2  and  N  =  4. 
The  solid  and  dashed  lines  correspond  to  T^-  =  Tc,  and  Tj<  =  10TC,  respec- 
tively. The  concentration  of  impurities  in  all  cases  is  n  =  0.4. 
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Figure  3.4  Conduction  electron  density  of  states  in  the  Coqblin-Schrieffer 
limit  for  N  =  2  and  several  values  of  hybridization  T;  T#  =  T  in  the  CS  limit. 
The  impurity  concentration  is  n  =  0.04. 
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Figure  3.5  The  position  and  spectral  weight  of  the  bound  states  in  the  gap 
for  an  s-wave  superconductor  with  Kondo  impurities.  The  solid  (dashed)  line 
is  the  location  (spectral  weight)  of  bound  states.  Only  one  of  the  bound  states 
is  indicated  here,  the  other  one  is  located  at  positive  energies,  symmetrically 
with  respect  to  the  gap  center.  Circles  (triangles)  refer  to  the  position  (spectral 
weight)  obtained  from  an  NRG  calculation  by  Shiba  et  al.  [52] 


Monte  Carlo  study  by  Jarrell  et  al.    [  63]  confirms  the  overall  dependence  of 
ujb  on  TK/TC. 

Bound  states  correspond  to  the  poles  of  the  T-matrix  for  conduction  elec- 
trons, T(uj)  =  V  Gt(u))  and  are  given  by 

1/2 


\uJ\ 


tk2  + 


rV 


T\u\ 


(3.15) 


We  compare  the  solution  of  equation  (3.15)  with  the  NRG  result  in  Figure 
3.5  .  There  is  a  good  agreement  for  Tr  >  Tc.  For  Tjj  ^>  Tc  the  position  of 
bound  states  ug/A  is  a  quadratic  function  of  Tx/Tc,  |o;^/A|  ~  1  —  2A  /T  %, 
and  agrees  with  NRG.  The  discrepancy  between  the  NRG  result  and  our  cal- 
culation for  Tfc-  <  Tc  is  not  surprising  since  the  slave-boson  theory  fails  in 
the  high  temperature  regime.  The  slave-boson  mean-field  amplitude  vanishes 
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around  T  ~  Tr  and  the  theory  is  unable  to  describe  the  crossover  to  temper- 
atures above  Tr.  In  Figure  3.5  we  also  show  the  spectral  weight  of  the  bound 
states  in  the  gap  which  again  agree  with  the  NRG  calculation  for  Tr  3>  Tc 

[58]. 

3.2.3  Specific  Heat  Jump 

To  calculate  the  specific  heat  jump  we  use  equations  (3.4),  (3.6)  and  (3.10), 
(3.11),  and  expand  the  gap  equation  near  Tc  in  terms  of  A/T, 


ln^  =  -^An  +  l/2)(l\a/BM)+^^2  +  lh{Sr)   '  ^ 


and  E(ujn)  =  1  +  a/B(u>n).  After  expanding  the  free  energy  to  fourth  order  in 
A  we  obtain 


^c)-C„W  =  ^-g9g^|.  (3.18) 

As  can  be  seen  in  Figure  3.6  ,  the  dependence  of  AC/ACq  on  rc/Tco,  for 

Tj{<Tcg,  is  qualitatively  different  from  the  Fermi  liquid  regime.  In  that  limit, 
C*  ~  1  —  l/ln(Tft-/2Tco),  see  Figure  3.7  .  Ichinose  [55]  and  Sakurai  [  64]  have 
obtained  a  somewhat  different  result,  C*  ~  1  —  1/ln  {Tj^/Tcq). 
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Figure     3.6     Specific  heat  jump  as  a  function  of  Tc/Tcq.   The  result  of  the 
Abrikosov-Gor'kov  theory  is  included  for  comparison. 
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Figure  3.7  The  derivative  of  the  specific  heat  jump  evaluated  at  zero  im- 
purity concentration  as  a  function  of  Tj^/Tc  for  an  s-wave  superconductor. 
The  dash-dotted  line  shows  the  asymptotic  behavior,  1  —  1/\ii(Tj^/2Tcq).  The 
dashed  line  is  the  asymptotic  form  at  T^/Tc  — >  oo  found  in  Refs.  [49]  and  [58], 

C*  ~  1  -  l/ln2(7V/Tco). 
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3.2.4  Penetration  Depth 

If  an  electromagnetic  wave  of  frequency  u  is  normally  incident  on  plane  su- 
perconducting     surface,       the      current      response      may      be      written 
j(q,  £2)  =  -K(<\,  £2)A(q,  £2),  where  A  is  the  applied  vector  potential.  The  pen- 


etration depth,  A 


4tt 


A'(0,0) 


tion  of  Skalski  et  al.  [  65], 


-1/2 


,  can  be  calculated  following  the  deriva- 


A 


mc 


_A_y/2 


n>0  ((1 


((1  +  (an/A)2)l/2  +  Qr/A^)(l  +  (w„/A)2) 


-1/2 


(3.19) 


where  A  =  lj2  +  T^  +  2ujnTCon/ \J A2  +  u2.  The  zero  temperature  penetration 
depth  is  given  in  the  strong  coupling  limit  by 


A(0)  a 

~  1  + 


l  +  nL)^NTco?2 


(3.20) 


AL(0)  "        '   2T2.  V^   '   4A7  "  16A0T2,' 

which  for  N  =  2  is  equal  to  1.071  in  the  Coqblin-Schrieffer  limit.  A  numerical 
calculation  of  the  penetration  depth  shows  that  the  derivative  d(X{0)/\L(0))/dfi, 

evaluated  at  n  =  0,  is  a  weakly  varying  function  of  Tk/Tcq  for  1  <  Tk/Tcq  < 
oo.  The  penetration  depth  is  much  less  affected  by  impurity  doping  on  the  high 
temperature  side,  where  d(\(Q)/\L(0))/dn\n=0  reaches  values  much  smaller 
than  those  in  the  strong  coupling  regime.  We  now  see  that  the  superfluid 
density  at  T  =  0  does  not  scale  with  the  Tc  suppression  in  the  Fermi  liquid 
regime.  While  the  slope  of  the  initial  Tc  suppression  goes  to  0  in  the  limit 
Tk/TcQ  -»  oo,  the  slope  of  the  initial  increase  of  the  T  =  0  penetration  depth 
remains  large. 


38 

3.3  Unconventional  Superconductors 

The  problem  of  an  Anderson  impurity  in  an  unconventional  superconductor 
is  of  considerable  interest  for  several  reasons.  First,  it  serves  as  a  testing 
ground  for  the  ideas  of  Withoff  and  Fradkin,  who  argued  that  in  the  analogous 
Kondo  problem  with  power  law  conduction  electron  density  of  states,  N(u>)  = 
C|u;|r,  there  existed  a  critical  coupling  below  which  impurities  are  effectively 
decoupled  from  the  conduction  band  [3,  66].  We  showed  in  the  preceding 
chapter  that  in  an  unconventional  superconductor  a  critical  coupling  indeed 
exists  for  the  1-impurity  problem,  but  that  for  a  finite  density  of  impurities 
there  is  always  a  finite  density  of  conduction  electron  states  at  the  Fermi  level, 
provided  N  =  2  or  if  r  <  1.  This  is  consistent  with  phenomenological  studies 
of  potential  scattering  in  unconventional  superconductors  [5,6].  This  does  not 
exclude  the  possibility  of  observing  some  vestige  of  this  transition  as  the  Kondo 
temperature  becomes  quite  small,  however.  Secondly,  Kondo  impurities  in 
unconventional  superconductors  have  been  proposed  as  analogues  of  defects  in 
Kondo  lattices  [  67],  with  the  argument  that  a  vacancy  in  a  Kondo  lattice  may 
induce  a  relative  phase  shift  close  to  n/2. 

We  chose  to  do  our  calculations  for  the  axial,  A(k)  =  A.Q(kx  +  iky),  and 
polar,  A(k)  =  Aq^z,  states  for  simplicity.  These  two  states,  nominally  p-wave 
pairing  states  over  a  spherical  Fermi  surface  in  3D,  are  quite  generally  rep- 
resentative of  two  classes  of  order  parameters.  States  with  order  parameters 
vanishing  at  points  on  the  Fermi  surface,  like  the  axial  state,  have  low  temper- 
ature properties  associated  with  a  quasiparticle  density  of  states  N(u>)  ~  u;2, 
whereas  the  states  with  lines  of  nodes,  like  the  polar  state,  correspond  to 
N(u)  ~  u.  More  complicated  order  parameters,  such  as  those  having  a  d-wave 
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symmetry,  will  have  similar  properties  at  low  temperatures,  since  the  main 
factor  determining  the  low-T  properties  is  the  order  parameter  topology,  i.e. 
whether  there  are  points  or  lines  of  zeros  of  the  order  parameter.  In  the  case 
of  points  (lines)  of  nodes,  the  low-temperature  specific  heat  of  pure  supercon- 
ductors is  proportional  to  T  (T).  The  deviation  of  the  penetration  depth  AA 
from  its  zero  temperature  value  A(0)  along  the  main  axes  of  symmetry  of  the 
order  parameter  is  either  ~  T  or  ~  T  in  the  axial  state,  according  to  the 
direction  of  current  flow  [  68],  whereas  for  the  polar  state  it  is  either  ~  T  or 
~  T3.  The  presence  of  strong  impurity  scattering  complicates  the  picture  and 
these  power  laws  do  not  hold  in  general. 

3.3.1  Suppression  of  the  Critical  Temperature 

Due  to  the  absence  of  the  Anderson  theorem  for  p-wave-like  or  d-wave-like 
superconductors,  the  influence  of  impurity  scattering  on  the  critical  tempera- 
ture is  qualitatively  different  from  that  of  s- wave-like  superconductors  [4,  69]. 
For  the  unconventional  states  of  interest,  there  are  no  off-diagonal  corrections 
to  the  superconducting  Green's  function,  A  =  A. 

In  this  case  Tc  for  either  type  of  state  is  determined  by 

The  initial  Tc-suppression  is  then  given  by 

(±dlc\  N   y,        £K  +  T) 

\TcO  dn)n=0         4tt2  Zj  (n  +  l/2)2B(cony 

and  approaches  -NT2/8TJ-,  in  the  limit  T^/TcQ  ->■  oo,  see  Figure  3.8  .  Note 
that  CJn/A  — >  aT/ATps  ,  as  Tc  — >  0  ,  and  as  a  consequence  the  critical  concen- 
tration nc  is  finite. 
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Figure  3.8  Critical  temperature  vs.  impurity  concentration  for  unconven- 
tional states  considered  in  this  work.  The  initial  slope  at  Tcq  is  shown  in  the 
inset. 


3.3.2  Density  of  States 

The  effect  of  Kondo  impurities  on  the  density  of  states  is  shown  in  Figures 
3.9  and  3.10  .  As  in  the  s-wave  case,  the  resonant  states  move  towards  the  edges 
of  the  gap  as  Tj^/Tc  increases,  except  in  the  special  case  N  =  2,  where  they 
are  pinned  at  the  gap  center.  In  all  cases  the  impurity  bands  are  broadened 
relative  to  the  s-wave  case  by  the  continuum  in  which  they  are  embedded. 
In  the  limit  Tr  — >  oo,  the  current  theory  coincides  with  the  results  given  by 
phenomenological  T-matrix  treatments  [5,6]  for  N  =  2. 


3.3.3  Specific  Heat 

To  obtain  the  specific  heat  we  differentiate  the  entropy  with  respect  to 
temperature  Cv  =  TdS/dT.  The  entropy  of  conduction  electron  quasiparticles 
renormalized  by  impurities  is  given  by 
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Figure  3.9  Conduction  electron  and  impurity  spectral  functions  in  the 
Kondo  limit  in  a  polar  state  for  N  =  2  and  N  =  4.  The  solid  and  dashed 
lines  correspond  to  Tr  =  Tc,  and  T^-  =  10TC,  respectively.  The  concentration 
of  impurities  is  n  =  0.4. 


roc 

S  =  -kB  /      duN(u)[f\nf  +  (1  -  /)ln(l  -  /)], 
JO 


(3.23) 


where  /  =  f(u)  is  the  Fermi  function.  Note  that  the  density  of  states  N(u)  is 
calculated  self-consistently,  using  equations  (3.2),  (3.3),  and  (2.20). 

The  presence  of  resonances  at  low  energies  leads  to  pronounced  features 
in  the  low-T  specific  heat,  as  shown  in  Figure  3.11  .  For  Tr  3>  Tc,  N  =  2, 
these  features  are  identical  to  those  predicted  by  the  phenomenological  theory 
of  Refs.  [5,6,  70].  For  smaller  Tj(/Tc,  the  resonance  sharpens,  as  is  evident 
from,  e.g.,  Figure  3.10  .  In  this  case  resonances  in  the  low-temperature  specific 
heat  may  be  quite  dramatic.  Figure  3.12  shows  the  temperature  dependence  of 
C/T  forN  =  2  in  the  axial  state  for  the  case  when  the  resonance  is  very  narrow, 
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Figure  3.10  Density  of  states  for  conduction  electrons  in  axial  (at  top)  and 
polar  (at  bottom)  state  for  N  =  2  and  T^/Tcq  =  0.3  (solid  line),  1  (dashed 
line),  and  20  (dash-dotted  line).  The  inset  shows  the  impurity  spectrum  for 
Tk/TqQ  =  20  in  both  states.  The  scaled  impurity  concentration  is  n  =  0.2. 


r  <C  Tc  and  just  above  the  Fermi  surface,  e  t  <C  Tc.  Such  pronounced  features 
are  possible  only  when  the  bare  impurity  level  is  close  to  the  Fermi  surface, 
and  the  hybridization  is  weak.  They  will  be  sharper  for  superconducting  states 
with  larger  exponent  r  in  the  unperturbed  density  of  states  N(u>)  ~  C|o;|r  at 
low  u.  These  anomalies  may  be  observed  experimentally  at  sufficiently  low 
temperatures.  In  this  context,  it  is  interesting  to  note  that  a  sharp  peak  in 
C/T  has  been  observed  in  the  heavy  fermion  superconductor  UPt%  at  18  raK 
[  71].  This  peak  is  present  at  roughly  the  same  position  also  in  the  normal 
state  at  magnetic  fields  B  >  Bc2,  as  might  be  expected  in  a  situation  where 
the  Kondo  temperature  is  significantly  smaller  than  the  critical  temperature. 
If  such  an  interpretation  of  the  measurement  of  Schuberth  et  al.  in  these  terms 
is  correct,  we  would  expect  the  size  of  the  peak  to  scale  with  other  measures 
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Figure  3.11  The  low-T  part  of  C/T  in  the  polar  state  for  N  =  2  in  the 
Kondo  limit.  The  inset  shows  C/T  for  n  =  0.01nc  over  the  full  temperature 
range. 


of  the  defect  concentration,  such  as  Tc  the  size  of  the  specific  heat  jump,  in 
different  samples.  We  note,  however,  that  fields  of  order  1  Tesla  would  normally 
destroy  a  many-body  resonance  of  the  usual  magnetic  type. 

3.3.4  Specific  Heat  Jump 

The  specific  heat  jump  at  Tc  can  be  found  by  the  method  already  mentioned 
in  the  preceding  section, 


CS(TC)  -  Cn(Tc) 


8n2N0Tc 


x  _  y^  ( G(un) 

^0i((n  +  l/2)%)+a/%))2 


(3.24) 
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Figure     3.12     The  low-T  part  of  C /T  in  the  axial  state.    Note  the  sharp 
resonance  in  the  density  of  states  at  low  energies  (the  inset). 


where 


h-Y. 


b  +  cH{un) 


(3.25) 


n>Q[(n  +  l/2)E(uJn)  +  a/B(u;n)}6 
with  b  =  4/5  (3/5),  and  c  =  2/3  (1/3),  for  the  axial  (polar)  state,  and 
H{un)  =  aT  [2(un  +  T)2/B(un)  -  l]  /B(un)  ,  G{un)  =  2au:n{LOn  +  T){n  + 
1/2  +  T/2irTc)/B2(u;n)  +  aT/2nTcB(un)  . 

The  derivative  of  the  specific  heat  jump,  C*,  is  shown  in  Figure  3.13  . 


3.3.5  Penetration  Depth 

We  now  discuss  the  low-temperature  response  functions  of  the  supercon- 
ductor in  the  presence  of  Kondo  impurities,  which  have  not  to  our  knowledge 
been  previously  calculated  in  the  strongly  interacting  Fermi  liquid  regime  of 
interest.   The  effect  of  Kondo  impurities  on  the  electromagnetic  response  has 


45 


5.0 


4.0 


6  3.0 


2.0 


1.0 


■2.0      -1.0        0.0         1.0 
logi0TK/Tc 


2.0 


Figure  3.13  The  derivative  of  the  specific  heat  jump  evaluated  at  zero 
impurity  concentration  as  a  function  of  T^/Tc  for  the  axial  and  polar  super- 
conducting states. 


been  calculated  in  the  phenomenological  T-matrix  approach  mentioned  above, 
and  used  to  analyze  experiments  on  heavy  fermion  and  high-Tc  supercon- 
ductors with  impurities,  but  no  microscopic  theory  is  available.   The  London 


4n 


A'(0,0) 


1-1/2 


,  where  A'(0,  0)  is 


penetration  depth  is  obtained  from  A  = 

the  electromagnetic  response  kernel.  The  kernel  is  given  by  the  linear  response 

formula 


3f>2y    roc         a  a 

KV(0Mm)  =  -- /      deY^r^{kikjG(k1Un)G(k1un-Qm)yk.     (3.26) 
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In  the  static  limit,  flm  — »•  0  ,  the  kernel  becomes 
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Figure    3.14     The  inverse  square  of  the  penetration  depth  vs.  temperature 
at  several  concentrations  of  impurities  in  the  polar  state. 


where  the  integral  represents  the  angular  average. 


Here  we  specialize  again  to  the  N  =  2  case.  The  temperature  dependence 
of  (A  )  in  the  polar  state,  which  corresponds  to  the  component  of  the  su- 
perfluid  density  within  the  plane  containing  the  line  of  zeroes  of  the  order 
parameter,  is  shown  in  Figure  3.14  .  The  low-T  behavior  changes  from  linear 
to  quadratic  upon  doping.  A  similar  result  was  obtained  earlier  [68,  72]  within 
the  phenomenological  theories.  The  results  for  T^  ~  Tc  are  qualitatively  very 
similar  to  those  shown.  Finally,  A  at  T  =  0  along  the  main  axes  of  symmetry 
for  the  axial  and  the  polar  state  as  a  function  of  impurity  concentration  is  pre- 
sented in  Figure  3.15  .  We  note  that  the  largest  component  of  the  penetration 
depth  scales  in  the  case  of  line  nodes  as  A(n)/A^(0)  —  1  ~  n1'2  at  low  con- 
centrations [  73],  and  A     (n)  ~  log(nc/n)  close  to  the  critical  concentration. 
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Figure  3.15  The  inverse  square  of  the  T  =  0  penetration  depth  in  the  two 
principal  directions  as  a  function  of  impurity  concentration  for  the  polar  (full 
lines)  and  the  axial  state  (dashed  lines). 


Similar  suppression  of  the  superfluid  density  was  calculated  for  a  d-wave  order 
parameter  within  a  phenomenological  model  in  Ref.  [  74]  .  It  is  clear  that 
both  the  concentration  and  temperature  dependence  of  the  penetration  depth 
components  may  be  important  tests  of  gap  anisotropy. 


CHAPTER  4 

ANISOTROPIC  CONVENTIONAL 

SUPERCONDUCTORS  VS. 

UNCONVENTIONAL  SUPERCONDUCTORS 


4.1  Experimental  Motivation 
Recent  high  resolution  angle-resolved  photoemission  (ARPES)  experiments 
on  Bi-2212  have  been  interpreted  in  terms  of  a  highly  anisotropic  order  param- 
eter with  a  large  gap  in  the  (ir,  0)  direction  and  two  lines  of  nodes  near  a  small 
gap  in  the  (7r,  tt)  direction  [  75,  76]. 


a)  b)  c) 

Figure  4.1  Order  parameters  plotted  over  tetragonal  Fermi  surface:  a) 
dx2_y2  state,  b)  extended  s-wave  state;  c)  s-wave  state  with  gap  minima  at 
Fermi  surface. 


While  alternative  explanations  consistent  with  an  order  parameter  A^ 
which  reduces  the  symmetry  of  the  Fermi  surface  (e.g.  d-wave,  see  Figure  4.1 
a)  have  been  put  forward  [  77],  it  is  interesting  to  consider  the  consequences 
of  the  simpler  suggestion  that  Aj.  has  the  full  symmetry  of  the  crystal,  but 
changes  sign  (see  Figure  4.1  b).  Fehrenbacher  and  Norman  [  78]  have  consid- 
ered the  effects  of  potential  scattering  on  such  a  state,  following  earlier  work 
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on  anisotropic-s  states  with  nodes  [  79],  and  shown  that  small  concentrations 
of  impurities  can  lead  to  "gapless"  behavior  in  the  density  of  states  (i.e.,  with 
residual  density  of  states  N(0)  >  0),  followed  by  the  opening  of  an  actual 
induced  gap  in  the  quasiparticle  spectrum  as  the  concentration  is  increased 
further. 

The  prediction  of  a  range  of  impurity  concentrations  over  which  "gapless" 
behavior  is  predicted  is  important  because  microwave  and  NMR  experiments, 
particularly  on  Zn  and  Ni-doped  YBCO  crystals  [  80,  81,  82,  83]  have  shown 
evidence  of  low-temperature  thermodynamic  properties  reflecting  the  existence 
of  a  residual  density  of  states.  They  have  been  interpreted  most  often  in  terms 
of  a  d-wave  pairing  scenario  [73,  84,  85,  86,  87],  where  the  well-known  effect 
of  dirt  is  to  induce  a  residual  density  of  states  for  infinitesimal  concentrations 
[69].  The  intriguing  possibility  raised  by  the  ARPES  experiment  [75,76]  is  that 
many  of  the  observations  of  "gapless"  behavior  can  be  equally  well  explained 
by  an  anisotropic  state  with  extended-s  symmetry.  There  are  some  difficulties 
associated  with  this  interpretation  even  if  it  can  successfully  account  for  the 
microwave  and  NMR  data.  The  most  important  of  these  is  the  set  of  SQUID 
experiments  [12,13,14,15]  indicating  that  the  order  parameter  Aj.  changes  sign 
under  a  ir/2  rotation.  While  the  extended-s  order  parameter  discussed  here 
indeed  changes  sign  over  the  Fermi  surface,  the  symmetry  of  the  expected 
tunneling  currents  is  not,  within  the  simplest  theory,  consistent  with  the  ob- 
servations reported  [  88,  89].  Furthermore,  the  existence  of  a  small  gap  in  the 
(tt,  tv)  direction  in  BSCCO  has  been  questioned  by  at  least  one  other  photoe- 
mission  group  claiming  similar  angular  and  energy  resolution  [  90].  We  do  not 
address  any  of  these  discrepancies  in  this  work. 
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Figure  4.2  ARPES  determination  of  BSCCO  energy  gap  as  function  of  angle 
0  around  Fermi  surface.  Data  from  Ref.  1.  Solid  line:  A^  =  Aq(|  cos2^|  —tjq), 
with  ??o  =  0.25. 


4.2  The  Model 

In  the  absence  of  definitive  answers  to  the  above  questions,  we  assume 
the  plausibility  of  the  Argonne-U.  Illinois  argument  and  investigate  the  conse- 
quences of  assuming  the  order  parameter  symmetry  identified  in  Refs.  [75,76] 
within  a  generalized  BCS  model. 

While  these  authors  pointed  out  that  their  data  were  consistent  with  an 
order  parameter  over  a  realistic  BSCCO  Fermi  surface  with  symmetry  Aj.  ~ 
cos  kx  cos  ky ,  we  work  here  with  an  even  simpler  model  with  cylindrical  Fermi 
surface,  in  which  case  a  rather  good  fit  to  the  data  can  be  obtained  by  assuming 
Afc  =  Aq(|cos2</>|  —  ?jq),  with  Ao=35meV  and  i]q  =  0.25,  as  shown  in  Figure 
4.2  .  With  this  set  of  parameters,  a  BCS  weak-coupling  approach  yields  the 
gap  magnitude  ratio  Aq/Tc  =  2.92,  see  Figure  4.3  . 
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Figure    4.3     The  ratio  Aq/Tco  as  a  function  of  tjq. 


Note  that  770  controls  the  angular  range  over  which  the  s-wave  order  param- 
eter is  negative.  The  case  tjq  =  0  corresponds  to  a  somewhat  pathological  state 
of  type  c),  while  770  =  2/ir  corresponds  to  equal  positive  and  negative  weights, 
(Afc)  =  0'  where  (...)  is  a  Fermi  surface  average.  When  7/0  =  0,  the  energy  of 
quasiparticle  excitations  is  identical  to  that  of  the  d-wave  state  dx2_yi.  Other 
choices  of  basis  functions  may  result  in  a  somewhat  better  overall  fit  and  avoid 
the  cusp  at  0  =  7r/4,  but  we  do  not  expect  these  details  to  affect  our  qualitative 
conclusions. 

Measurements  of  penetration  depth,  angle-resolved  photoemission,  thermal 
conductivity,  and  nuclear  magnetic  relaxation  experiments  provide  evidence  for 
gap  nodes,  but  do  not  determine  if  the  order  parameter  changes  sign.  This  is 
because  experiments  of  this  kind  are  normally  assumed  to  measure  properties 
sensitive  to  the  order  parameter  Ak  only  through  the  Bogoliubov  quasiparticle 


spectrum,  Ek  =  Jg*  +  \Ak\2.    Since  Ek  does  not  depend  on  the  sign  of  the 
order  parameter,  all  three  states  shown  in  Figure  4.1  have  similar  (linear) 
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Figure  4.4     Density  of  states  for  the  superconducting  state  with  an  extended- 
s  order  parameter,  Aj.  =  Aq(|  cos20|  —  j]q),  with  770  =  0.25. 


low-energy  density  of  single-particle  states  N(u),  and  all  properties  deriving 
directly  therefrom. 

Figure  4.4  shows  density  of  states  for  the  state  of  the  type  4.1  c.  Note  the 
low-energy  gap  feature  corresponding  to  the  smaller  gap  maximum. 

Systematic  impurity-doping  studies  of  transport  properties  can  in  principle 
help  to  identify  the  order  parameter  symmetry,  however.  Unconventional  states 
like  the  dx2_y2  state  (Figure  4.1  a)  possess  nodes  on  the  Fermi  surface  for 
symmetry  reasons,  whereas  their  existence  is  "accidental"  in  the  s-wave  case 
(Figure  4.1  b-c).  In  the  d-wave  case,  supression  of  gap  variation  in  k-space  can 
only  result  in  an  overall  supression  of  the  order  parameter  magnitude,  whereas 
in  s-wave  cases  b)  and  c),  it  results  eventually  in  the  elimination  of  the  nodes. 
In  many-body  language,  the  crucial  point  is  that  the  off-diagonal  impurity  self- 
energy  Si  is  nonzero  in  s-wave  states  4.1  a  and  4.1  b,  but  vanishes  in  the  d- 
wave  state  4.1  c.  In  the  isotropic  s-wave  state,  a  large  Si  cancels  the  diagonal 
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self-energy  So  [1],  leading  to  Anderson's  theorem  [2],  whereas  in  extended-s 
states  this  cancellation  is  only  partial.  A  large  Si  furthermore  prevents  the 
formation  of  a  scattering  resonance  at  the  Fermi  surface,  leading  to  clean  limit 
low-frequency  transport  coefficients  which  in  extended-s  states  for  T  -C  Tc 
vary  relatively  weakly  with  temperature.  By  contrast,  in  the  d-wave  (or  other 
unconventional  state  where  Si  =  0),  resonant  scattering  leads  to  transport 
properties  which  vary  strongly  in  temperature  for  T  <C  Tc  [  91,5,6].  In  this 
sense  the  characteristics  attributed  to  extended-s  states  here  are  also  valid 
qualitatively  for  any  mechanism  which  generates  an  off-diagonal  self-energy  in 
a  d-wave  state  as  well,  e.g.  tetragonal  symmetry  breaking  or  locally  induced 
s-wave  components  due  to  impurity  potentials. 

The  disorder-averaged  matrix  propagator  describing  any  of  the  states  a)-c) 
above  is  written 

9{k,un)  =  — — ~s — -= — —  (4.1) 

^2  -  el- 1  Ak  i2 

where  the  t1  are  the  Pauli  matrices  and  Aj.  is  a  renormalized  unitary  order 
parameter  in  particle-hole  and  spin  space.  The  renormalized  quantities  are 
given  by  u  =  u  -  Eo(u>),  4  =  4  +  -3M>  and  Ak  =  Ak  +  Ei(o;),  where  the 
self-energy  due  to  s-wave  impurity  scattering  has  been  expanded  S  =  SjT1. 
The  relevant  self-energies  are  given  in  a  self-consistent  /-matrix  approximation 
[5,6]  by 

^°        c2  +  Gl2  _  Gq2  '   -1  -  c2  +  Gl2  _  Gq2  •  W 

where  T  =  nin/^No)  is  a  scattering  rate  depending  only  on  the  concentration 
of  defects  ??,?;,  the  electron  density  n,  and  the  density  of  states  at  the  Fermi 
level,  Nq,  and  we  have  defined  Ga  =  (i/2TrNQ)T,i.Tr[Tag].  The  strength  of  an 
individual  scattering  event  is  characterized  by  the  cotangent  of  the  scattering 
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phase  shift,  c.   The  Born  limit  corresponds  to  c  >  1,  so  that  So  ~  ±Tj^Go, 

i  i 

while  the  unitarity  limit  corresponds  to  c  =  0.  We  have  defined  the  normal- 
state  impurity  scattering  rate  as  Tj\f  =  T/(l  +  c  );  note  that  in  the  high-Tc 
cuprates  the  total  scattering  rate  at  Tc  includes  inelastic  scattering  and  is 
expected  to  be  much  larger  for  clean  samples. 

4.3  Density  of  States 

A  crucial  feature  of  the  physics  of  d-wave  superconductors  is  that  an 
infinitesimal  concentration  of  impurities  produces  a  finite  density  of  states 
N(0)  >  0  at  the  Fermi  level  [69],  leading  to  temperature  dependences  character- 
istic of  the  normal  state  in  all  transport  quantities.  Solving  the  self-consistency 
equations  at  u  =  0  for  the  extended-s  wave  state  under  consideration  leads  im- 
mediately to  the  conclusion  that  such  "gap less"  behavior  is  possible  only  for  a 
range  of  scattering  rates  T  <  Tc,  however  [78]. 

This  is  illustrated  in  Figure  4.5  .  A  low  frequency  expansion  in  the  gapless 
regime  yields  Tc/Tc0  ~  %(1  +  c2). 

Can  the  same  experiments  which  seem  to  fit  the  "dirty  d-wave"  scenario 
also  be  explained  by  extended-s  states?  The  difficulty  is  how  to  fix  the  actual 
impurity  scattering  rate,  T,  given  the  known  concentration.  One  way  is  to 
attribute  the  additional  extrapolated  T  — >  0  resistivity  to  impurity  scattering, 
such  that  the  elastic  and  inelastic  rates  add  incoherently.  If  one  attempts  such 
an  analysis  for  Zn-doped  YBCO  crystals  using  results  from,  e.g.,  Chien  et  al. 
[  92],  one  finds  that  T/Tc  ~  0.3  -  0.5  per  1%  Zn.  Since  Zn  doping  studies  of 
YBCO  indicate  gapless  behavior  up  to  several  per  cent  Zn,  it  would  appear 
that  an  extended-s  picture  is  plausible  for  YBCO  only  if  one  assumes  770  close 
to  the  critical  value  2/n  for  which  (Aj.)  =  0.    We  are  not  aware  of  similar 
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Figure  4.5  Residual  density  of  states  N(Q)/Nq  in  extended-s  state  vs.  nor- 
malized scattering  rate  r/Tco  for  c  =  0.  Dashed  line:  ?y0  =  °-25  obtained  from 
fit  to  data  of  Ding  et  al. 


doping  studies  in  single  crystal  BSCCO,  but  assuming  for  the  moment  that  Zn 
scatters  equally  strongly  in  this  material,  we  see  from  the  dashed  line  in  Figure 
4.5  that  a  well-developed  gap  ^Iq  in  the  excitation  spectrum  should  be  induced 
in  BSCCO  by  a  few  per  cent  Zn  doping.  Note  that  the  results  shown  in  Figure 

4.5  are  not  sensitive  to  changes  in  the  scattering  phase  shift  8q.  In  the  special 
case,  tiq  =  0,  a  simple  estimate  shows  that  for  small  scattering  rates,  Q,q  ~  T 
(r^v  in  Born  limit). 

In  the  dirty  limit  T  -)■  oo,  the  s-wave  superconductor  becomes  isotropic 
with  a  BCS  density  of  states  N(uj)  =  Re  uj/Juj2  -  Aavg2,  as  shown  in  Figure 

4.6  .  In  contrast  to  a  d-wave  superconductor,  the  self-energies  obtained  in  the 
Born  approximation  and  in  the  resonant  scattering  limit  are  almost  equivalent 
in  the  highly  anisotropic  s-wave  system,  if  770  <  AaUi,/Ao.  This  insensitivity  to 
larger  phase  shifts  arises  because  of  off-diagonal  self-energy  corrections  which 
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Figure  4.6  Normalized  density  of  states  N(lj)/Nq  for  s-  and  d-wave  order 
parameters  vs.  reduced  frequency  w/Aq,  shown  for  various  potential  scattering 
rates  Tyy/Ao  in  the  Born  approximation;  ?/q  =  0. 


prevent  the  occurence  of  poles  in  the  t-matrix,  c2  -  Gq2  +  G\2  ~  (9(1)  for 

all  c<l.    Densities  of  states  for  both  types  of  states  in  the  limit  of  resonant 
scattering  are  shown  in  Figure  4.7  . 
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Figure  4.7  Normalized  density  of  states  N(lu)/Nq  for  s-  and  d-wave  order 
parameters  vs.  reduced  frequency  u/Aq,  shown  for  various  potential  scattering 
rates  T/Aq  in  the  unitarity  limit,  c  =  0;  770  =  0. 


4.4  Critical  Temperature 

We  first  solve  the  Dyson  equation  for  the  renormalized  propagator  (4.1) 

together  with  the  gap  equation.  The  order  parameter  A^.  is  related  as  usual  to 

the  off-diagonal  propagator  as  A(k)  =  TJ2n^2k'  ^fcfc'Tr(Ti/2)g(&', wfl),  where 

V^i  =  Vj  s®d  s(k)®d  s(h')  is  the  phenomenological  pair  interaction  assumed. 
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In  the  d-wave  case,  Ak  =  AQd$d(k),  with  $d  =  cos  20,  yielding  a  simple 
equation  for  the  gap  maximum  Aq^, 


ruJD  du  $j 

\d~'  ----  (   /       dwtanh^-  Re  d      =  )  ,  (4.3) 


'■/ 


0  2  J&  -  A,2 


where  A^  =  V^Nq,  and  (...)  represents  an  angular  average  over  the  cylindrical 
Fermi  surface.  In  the  s-wave  case,  on  the  other  hand,  it  is  convenient  to  put 
Afc  =  Aavg  +  Ajfc,  where  Aavg  =  (A^.)  is  the  gap  average  over  the  Fermi 
surface.  When  impurities  are  added  to  the  system,  it  is  easy  to  check  that 
Aj.  =  Aavg  +  A]^,  determined  by 

\-l-{rj ■■*,"-*;A*^\  (4.4) 

Note  this  is  effectively  an  equation  for  Aavg(uj)  since  the  angular  variation  A^ 
is  given.  When  i]q  =  0  the  initial  slope  of  Tc  suppression,  dTc/dTN  =  -\7r/4, 
where  \  =  [($s2)  -  ($s)2]/($s)2  is  1  -8/7r2  for  the  s-wave  and  1  for  the  d-wave 
state  considered.  In  the  d-wave  case  the  critical  temperature  continues  to  drop 
rapidly  to  zero  at  a  critical  concentration  of  n^  =  7r2./VoTco/2e7,  whereas  the 
decrease  becomes  more  gradual  as  the  gap  is  smeared  out  in  the  s-wave  case, 
finally  varying  [  93,  94,  95]  as  Tc  ~  Tc0[l  -  \  ln(1.154r^/7rTc0)]. 

The  initial  suppression  of  Tc  depends  also  on  the  magnitude  of  770  and  the 
magnitude  of  the  spin  exchange  scattering  rate  Ts.  For  small  770  and  rs  ^  F 
the  qualitative  form  of  Tc-suppression  at  small  and  moderate  concentrations 
will  be  very  similar  to  the  one  shown  above.  For  large  concentrations,  such 
that  Tc/Tcq  <C  1,  the  Tc-suppression  by  the  magnetic  component  in  the  scat- 
tering will  be  the  dominant  effect,  if  Fs  >  0.  However,  one  may  not  treat  the 
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impurities  as  independent  from  one  another  when  the  impurity  concentration 
is  large. 

4.5  London  Penetration  Depth 
The  opening  of  the  energy  gap  with  increasing  impurity  concentration  is 
an  indelible  signature  of  s-wave  superconductivity.  It  will  obviously  give  rise  to 
activated  behavior  for  T  <C  Qq  in  a  wide  range  of  thermodynamic  properties, 
of  which  we  have  chosen  to  discuss  only  one  for  purposes  of  illustration,  the 
temperature-dependent  magnetic  penetration  depth.  For  the  model  states  and 
Fermi  surface  under  consideration,  this  may  be  expressed  as 

MQ2=/^tanh^/^Re(  *j  I,  (4.5) 

\MT))        J  2    J   2k         \  (0,2  _  a|)3/2  J ' 

where  Aq  is  the  pure  London  result  at  T  =  0. 

The  penetration  depth  in  a  d-wave  superconductor  (Figure  4.8  b)  is  known 
to  vary  as  A(T)  ~  Ao  +  C2±  at  the  lowest  temperatures  [68,  96,  97,73],  over  a 
temperature  range  which  widens  with  increasing  impurity  concentration.  The 
coefficient  C2  decreases,  as  T_1  in  the  Born  limit  and  T-1'2  in  the  resonant 
scattering  case.  In  the  anisotropic  s-wave  case  with  770  >  0  there  is  a  gradual 
change  in  the  temperature  dependence  of  the  penetration  depth  at  low  tem- 
peratures. The  deviation  from  the  linear  behavior  of  \{T)  grows  weakly  with 
increasing  impurity  concentration,  see  Figure  4.9  .  As  mentioned  earlier  the 
energy  gap  opens  up  at  Tc/Tc0  ~  r/0(l  +  c2). 

For  770  "C  1  the  corresponding  activated  behavior  in  the  anisotropic  s-wave 
case  is  easy  to  distinguish  from  the  d-wave  case  when  plotted  against  (T/Tc)  as 
also  shown  in  Figure  4.8  .  The  important  experimentally  relevant  signature  is 
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Figure  4.8  Temperature  dependence  of  normalized  magnetic  penetration 
depth  (Ao/A(T))2  for  s-  and  d-wave  order  parameters  (?/q  =  0)vs.  reduced 
temperature  (T/Tc)2,  shown  for  various  potential  scattering  rates  T/Tco  in  the 
unitarity  limit,  c  =  0. 


of  course  not  simply  the  exponential  behavior,  but  the  increase  in  the  activation 
gap  with  impurity  concentration. 


4.6  Effect  of  Spin  Scattering 
A  simple  defect  like  a  vacancy  or  Zn  ion  in  the  CuC>2  plane  may  not  behave 
simply  as  a  potential  scatterer,  as  assumed  above.  In  the  presence  of  large  local 
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Figure  4.9  Temperature  dependence  of  normalized  magnetic  penetration 
depth  (Aq/A(T))  for  an  order  parameter  of  an  extended-s  type,  with  r/g  =  0.25, 
in  the  unitarity  limit,  c  =  0. 


Coulomb  interactions,  a  magnetic  moment  may  form  around  the  defect  site, 
giving  rise  to  spin-flip  scattering  of  conduction  electrons  [  98] .  This  poses  the 
most  serious  obstacle  for  the  direct  application  of  the  principle  distinguishing 
d-wave  from  anisotropic  s-wave  systems  outlined  above,  since  magnetic  scat- 
tering will  lead  to  gapless  superconductivity  as  in  the  usual  Abrikosov-Gor'kov 
theory.  Furthermore,  even  if  a  gap  remains,  strong  spin-flip  scattering  may 
lead  to  bound  states  within  it  [66,  99,  100,46]  which  may  give  rise  under  the 
proper  circumstances  to  a  residual  density  of  states  N(u>  — >  0)  as  in  the  d-wave 
case.  Here  we  investigate  the  competition  between  the  opening  of  the  energy 
gap  in  the  s-wave  state  due  to  potential  scattering  and  gapless  behavior  due 
to  magnetic  scattering.  To  this  end  we  add  a  term  JS  •  a  to  the  Hamiltonian, 
where  S  is  a  classical  spin  representing  the  impurity  and  a  is  the  conduction 
electron  spin  density,  and  study  the  system  in  an  average  t-matrix  approxima- 
tion analogous  to  the  one  applied  to  the  pure  potential  scattering  case.   The 
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Figure  4.10  Induced  energy  gap  normalized  to  clean  gap  maximum,  Q^/Aq, 
vs.  potential  scattering  rate  T^/Tcq  for  different  ratios,  TSN/Tj\i,  of  magnetic 
to  potential  scattering  rates,  tjq  =  0. 


self-energies  found  in  the  presence  of  both  types  of  scattering  reduce  in  the 
isotropic  s-wave  case  to  those  given  by  several  authors  [46,99,  101].  We  find 
that  until  the  dimensionless  exchange  JNq  becomes  of  0(1),  the  results  for  the 
s-wave  system  are  very  similar  to  those  obtained  in  the  simpler  Born  approx- 
imation, as  discussed  above,  provided  ??o  is  not  large.  The  self-energies  in  the 
Born  approximation  are  Eq  =  (T^  +  Tsn)Gq  and  £j  =  (— Tjy  +  ^%)G\  [1], 
where  TSN  =  n{J  S(S  +  1)ttNq.  For  ?7q  =  0,  the  induced  gap,  tig,  in  the  s-wave 
system  may  then  be  shown  to  vary  as  Qg  ~  T^  —  TSN  >  0,  but  the  effects  of 
self-consistency  rapidly  become  important  as  the  concentration  is  increased. 

In  Figure  4.10  ,  we  plot  fig  as  a  function  of  the  impurity  concentration 
through  the  parameter  T^r  for  various  assumptions  about  the  scattering  charac- 
ter of  the  impurity  ion,  where  the  quantity  T^/T^  specifies  the  relative  amount 
of  magnetic  scattering.  The  destruction  of  the  induced  gap  takes  place  because 
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the  system  becomes  insensitive  to  large  amounts  of  potential  scattering,  but 
magnetic  impurities  continue  to  break  pairs  even  at  large  concentrations.  The 
gap  is  nevertheless  found  to  persist  into  the  very  dirty  limit  even  for  systems 
where  the  magnetic  scattering  is  nearly  as  strong  as  the  potential  scattering. 

For  weak  spin  scattering,  the  bound  state  in  the  t-matrix  approximation 
is  found  to  lie  at  u  >  £Iq,  just  below  the  average  order  parameter  Aavg  deep 
in  the  continuum,  and  thus  plays  no  role.  Stronger  spin  scattering  does  not 
change  this  qualitative  behavior  at  low  concentrations  until  JNq  ~  1  when  the 
bound  state  lies  at  the  Fermi  level  in  the  classical  spin  approximation  [46].  In 
this  case  the  Kondo  effect,  neglected  here,  also  becomes  important.  It  is  known 
from  other  analyses  [  102]  that  the  bound  state  lies  near  the  Fermi  level,  and 
will  therefore  give  rise  to  a  residual  density  of  states  N(u>  ->  0),  only  when 
TK  —  Tc-  For  any  other  ratio  of  TK/TC,  the  bound  state  will  lie  at  an  energy 
corresponding  to  an  appreciable  fraction  of  the  average  gap  in  the  system,  and 
hence  be  irrelevant  for  our  purposes. 

Clearly  a  quantitative  estimate  of  the  relative  size  of  V  and  TSN  is  required 
to  decide  whether  spin  scattering  plays  a  role  in  real  high-Tc  materials  with 
simple  defects.  Walstedt  and  co-workers  estimated  JNq  ~  0.015  for  a  Zn  ion 
in  YBCO,  implying  that  Zn  is  a  nearly  pure  potential  scatterer  in  this  system 
[18].  On  the  other  hand,  Mahajan  et  al.  [19]  estimate  JNq  ~  0.45.  For  a  1%  Zn 
concentration,  a  magnetic  moment  of  0.36  \iB  for  Zn  in  fully  oxygenated  YBCO 
[19]  and  a  density  of  states  of  1.5/eV  [19],  we  find  T^  ~  1  x  10~4  eV.  From  the 
residual  resistivities  of  Zn-doped  YBCO  crystals  [92],  we  estimate  that  a  1%  Zn 
sample  corresponds  to  a  total  impurity  scattering  rate  of  T^p  ~  1  x  10~2ey, 
assuming  that  the  inelastic  and  elastic  contributions  to  the  scattering  rate 
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add  incoherently.  This  suggests  that  potential  scattering  must  dominate  the 
total  elastic  rate,  TSN  <C  I\  On  the  other  hand,  the  large  value  of  JNq  ~  0.45 
deduced  for  a  Zn  ion  [19]  would  mean  that  the  Kondo  effect  may  be  important, 
and  that  we  cannot  completely  rule  out  the  possibility  that  a  bound  state  sits 
very  close  to  the  Fermi  level.  Other  experimental  results  do  not  support  such 
high  values  for  the  magnetic  exchange  interaction  of  Zn  in  YBCO.  NMR  studies 
of  the  Osaka  group  [80,82]  led  to  the  conclusion  that  Ts  <C  T,  and  confirm  our 
estimates  of  Ts. 

4.7  Microwave  Conductivity 

Electrical  and  thermal  conductivities,  as  well  as  sound  attenuation,  will  be 
considerably  different  in  extended  s-  and  d-wave  states,  as  suggested  by  the 
following  argument.  Any  DC  transport  coefficient  L(T)  in  a  system  charac- 
terized by  well-defined  single-particle  excitations  will  vary  with  temperature 
roughly  as  L(T)  ~  JV(u>  ~  T)t(uj  ~  T),  where  N(u)  is  the  density  of  states 
and  t(lj)  the  relaxation  time.  In  the  clean  d-wave  case,  resonant  scattering 
gives  t  (u)  =  2^2q"(u)  ~  7V(u;)-1  up  to  logarithmic  corrections,  yielding 
L(T)  ~  T  .  In  the  Born  limit,  c  >1,  similar  arguments  yield  L  ~  const,  for 
d-wave  transport  coefficients. 

Impurity-limited  transport  in  the  extended  s-wave  state  will  be  qualita- 
tively similar  to  the  d-wave  case  if  the  scattering  is  weak,  c  ^>  1.  In  contrast 
to  the  d-wave  case,  however,  for  c  — >  0  (and  tjq  <gi  1)  resonant  scattering  does 

not  occur,  since  the  denominator  of  the  t-matrix,  c    —  Gq  +  G^  >  (1  —  ?/q) 
[79].  A  simple  low-T  estimate  accounting  for  nodal  quasiparticle  contributions 
gives  L  ~  Lq(1  —  2?7q)  as  T  -»  0.    The  exact  behavior  in  this  range  will  be 
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Figure  4.11  Conductivity  of  extended-s  state,  A0/Tc=2.92,  r/Tc=0.1,  c  = 
0,  770=0.1,0.25,  and  0.5  (solid  lines).  Also  shown,  T/Tc  =  0.01,  tjq  =  0.25 
(dashed  line). 


influenced  by  self-consistency  effects  and  the  leading  frequency  dependence  of 
the  t-matrix.  The  resultant  temperature  dependence  will  then  be  intermediate 
between  the  strong  and  weak  scattering  limits  of  d-wave  transport  coefficients. 

In  Refs.  [  103,  104],  expressions  for  the  complex  conductivity  of  an  anisotropic 
s-wave  superconductor  were  derived.  We  do  not  reproduce  these  rather  lengthy 
expressions  here,  but  merely  comment  that  a  fully  self-consistent  numerical 
evaluation  confirms  the  qualitative  picture  of  low-frequency  transport  described 
above. 

Some  typical  results  are  shown  in  Figure  4.11  for  a  model  in  which  in- 
elastic scattering  has  been  neglected  entirely.  The  limiting  conductivity  as 
T  ->  Tc  is  therefore  the  impurity  Drude  result,  a0  =  ne2/(2mTN),  which  for 
Tat  <  Tc  is  much  larger  than  the  actual  conductivity  in  the  cuprates  at  the 
transition,  a(Tc),  indicated  roughly  in  the  figure.  At  low  temperatures  T  <C  Tc, 
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inelastic  scattering  may  be  neglected  and  the  results  displayed  are  valid.  The 
most  important  qualitative  feature  of  the  results  is  that  for  tjq  <  2/tt,  the 
effective  limiting  value  of  the  conductivity  in  the  extended-s  state  is  nonzero 
and  generically  much  larger  than  cr(Tc)  in  clean  high-Tc  systems,  such  that 
a  =  a(T  — >  0)/cr(Tc)  3>  1.  This  residual  conductivity  diminishes  as  770  — »  2/n, 
when  the  result  should  be  qualitatively  similar  to  the  d-wave  conductivity  with 
resonant  scattering  due  to  the  vanishing  of  the  off-diagonal  self-energy.  Note, 
however,  that  for  generic  values  of  770,  e.g.  that  apparently  appropriate  for 
BSCCO,  the  residual  conductivity  scales  inversely  with  the  impurity  scattering 
rate  T,  in  contrast  to  the  resonant  d-wave  case  where  the  residual  conductivity 
is  independent  of  F  to  leading  order. 

As  shown  in  Figure  4.11  ,  the  temperature  dependence  of  the  conductivity 
for  most  of  the  low  temperature  range  may  mimic  a  linear  behavior.  This  is 
intriguing,  given  the  linear- T  conductivity  observed  in  microwave  experiments 
on  clean  YBCO  crystals  [81],  but  the  large  residual  conductivity  predicted 
(unless  770  is  close  to  2/n)  would  seem  to  be  inconsistent  with  recent  studies 
indicating  that  a(T  — >•  0)  is  very  small  in  twin-free  samples.  We  are  not  aware 
of  similar  experiments  on  high-quality  BSCCO  single  crystals.  One  final  point 
of  interest  is  that  the  form  of  o~(T)  is  only  weakly  dependent  on  disorder, 
as  shown  in  the  figure;  of  course  the  overall  conductivity  scale  o~q  depends 
inversely  on  impurity  concentration. 

The  results  in  the  Born  approximation  are  qualitatively  similar  for  both  s- 
and  d-wave  states,  as  can  be  seen  from  comparison  of  Figure  4.12  and  Figure  2 
of  Ref.  [104].  Note,  however,  that  the  clean  limit  regime  in  the  d-wave  state 
is  a  much  more  narrow  range  of  scattering  rates  than  in  the  s-wave  states. 
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Figure    4.12     Conductivity  of  extended  s-state  for  several  values  of  external 
frequency  Q  in  the  Born  approximation,  Aq/Tc  =  2.92,r/Tc0  =  0.1,  r/0  =  0.25. 
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Figure  4.13  Comparison  of  conductivity  of  clean  extended  s-wave  state 
(dashed  line),  t?0  =  0.25,  A0/Tc  =  2.92,  T/Tc  =  0.1,  c  =  0  with  dxi_  2  state, 
A/Tc  =  2.14,  r/Tc  =  0.1,  c  =  10, 1,  and  0. 
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To  illustrate  the  comparison  of  these  results  with  those  expected  for  a 
d-wave  superconductor,  we  plot  in  Figure  4.13  the  conductivity  in  a  dx2_y2 
state  in  the  resonant  (c  =  0)  and  Born  (c  =  10)  limits,  together  with  an 
example  of  intermediate  strength  scattering  (c  =  1)  chosen  to  give  the  same 
effective  a(T  — >  0)  as  the  extended-s  conductivity  for  770  =  0.25.  Even  in 
the  last  case  the  qualitative  differences  between  the  s-  and  d-wave  results  are 
manifest.  The  dependence  of  the  s-wave  result  on  the  scattering  phase  shift 
and  impurity  concentration  are  found  to  be  quite  weak,  except  in  the  case 
c  3>  1,  for  which  the  result  is  qualitatively  similar  to  the  d-wave  Born  result, 
since  the  denominator  of  the  t-matrix  becomes  irrelevant. 

Given  the  uncertainty  regarding  the  origin  of  the  residual  conductivity  in 
the  high-Tc  materials  it  is  perhaps  useful  to  give  estimates  for  other  transport 
coefficients.  For  example,  since  the  coherence  factors  are  essentially  the  same 
for  the  electronic  thermal  conductivity  nei/T  as  for  cr(T),  it  is  straightforward 
to  see  that  the  large  residual  conductivity  in  the  extended  s-wave  state  implies 
a  large  linear-T  term  in  the  thermal  conductivity  of  order  ^(1  —  27]q)^v  pT/Tj^f, 
where  vp  is  the  Fermi  velocity  and  7  is  the  normal  state  linear  specific  heat 
coefficient.  Such  a  term  has  not  been  observed  below  2A'  in  single  crystal 
BSCCO  [  105]  but  there  are  data  at  somewhat  higher  temperatures  which  are 
consistent  with  &q  =  nei(T  — )■  0)TC/ nei(Tc)T  ^>  1  [  106].  The  linear  term  in  k 
at  low  T  in  YBCO  appears  to  be  quite  small  [  107]. 

Very  recently  the  authors  of  Ref.  [  108]  have  pointed  out  that  the  "bump" 
in  BSCCO  gap  function  near  4>  =  n/4  may  be  a  superlattice  effect.  In  this  case 
available  experimental  data  on  BSCCO  seem  unlikely  to  be  compatible  with 
an  anisotropic  s-wave  state  of  the  type  considered  here  for  any  770- 


CHAPTER  5 

TWO-CHANNEL  KONDO  IMPURITIES 

IN  SUPERCONDUCTORS 

In  recent  years  models  of  fermion  systems  that  cannot  be  described  by  a 
Landau  Fermi  liquid  model  received  considerable  attention.  One  of  those  is  the 
two-channel  Kondo  impurity.  In  this  model  conduction  electrons  described  by 
the  flavor  a  and  spin  a  interact  with  impurity  spin  S  =  1/2.  It  was  pointed  out 
by  and  Nozieres  and  Blandin  [  109]  that  this  model  may  follow  from  a  realistic 
description  of  magnetic  impurities  in  a  metal  when  the  orbital  structure  of  the 
impurity  is  taken  into  account.  In  the  two-channel  model  the  impurity  spin  is 
overscreened  by  two  exactly  degenerate  channels  of  conduction  electrons.  The 
new  effective  spin  is  1/2.  This  spin  couples  antiferromagnetically  to  two  other 
electrons  from  the  two  channels,  and  this  process  continues  as  T  is  lowered.  The 
effective  coupling  of  the  conduction  band  grows  in  this  process  and  saturates 
at  large,  finite  values  at  an  intermediate  fixed  point.  The  NRG  approach 
shows  that  the  size  of  the  electron  cloud  around  the  impurity  diverges  as  1/T 
when  T  — >  0.  The  excitation  spectrum  of  this  system  is  not  a  Fermi  liquid. 
The  low  temperature  magnetic  susceptibility  and  the  specific  heat  coefficient, 
Jimp  =  Cimp/T,  are  logarithmically  divergent.  Other  properties  also  exhibit 
non-Fermi  liquid  behavior  below  the  crossover  temperature  Tg. 

The  experimental  realization  of  the  two-channel  Kondo  effect  in  uranium 
compounds  was  proposed  but  remains  controversial  [  110,  111,  112].  It  was 
argued  [110]  that  in  compounds  such  as  UBei3  and  Yi_xUxPd3  with  x  =  0.2, 
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the  U  ions  at  crystal  sites  with  cubic  symmetry  have  stable  5/  configuration 
and  J  =  4  total  angular  momentum.  The  crystal-field  split  J  =  4  multiplet 
would  have  a  nonmagnetic  I^  doublet  as  its  ground  state.  In  a  simplified 
model  this  doublet  is  coupled  to  an  excited  doublet  in  an  5/1  configuration, 
with  J  =  5/2  and  I^  symmetry,  via  hybridization  with  conduction  electron 
J  =  5/2,  Tg  partial  waves.  A  canonical  transformation  [31,110]  yields  an 
effective  interaction  between  the  T^  doublet  and  two  channels  of  conduction 
states  of  Fg  symmetry. 

Here  we  would  like  to  offer  another  test  of  the  multichannel  Kondo  effect. 
The  breaking  of  time-reversal  symmetry  involved  in  scattering  by  magnetic 
impurities  in  superconductors  has  important  consequences.  The  qualitative 
type  of  behavior  is  well-known  in  the  single  channel  problem  to  depend  on  the 
ratio  Tft-/Tcg,  where  Tj^  is  the  Kondo  temperature  and  Tcq  is  the  transition 
temperature  of  a  pure  superconductor.  The  divergence  of  magnetic  correlation 
length  and  the  existence  of  a  residual  magnetic  moment  as  T  — >  0  suggest  that 
in  the  multichannel  case  one  should  expect  a  very  different  type  of  interplay 
between  exchange  and  pairing  interactions. 

One  may  use  a  sensitivity  of  superconducting  correlations  to  multichannel 
exchange  interaction  as  an  additional  criterion  to  characterize  multichannel 
Kondo  behavior  and  distinguish  it  from  the  single  channel  case.  In  this  paper 
we  study  a  simplified  SU(N)xSU(k)  version  of  the  full  multichannel  problem 
in  the  NCA  approximation  [  113].  Here,  iV  is  the  orbital  degeneracy  of  the 
impurity  and  k  is  the  number  of  conduction  electron  channels  coupled  to  the 
impurity.  Although  this  approach  ignores  details  of  anisotropic  exchange  it 
remains  in  the  same  universality  class  [113].   In  particular  the  exponents  for 
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Figure  5.1  The  leading  order  diagram  for  the  anomalous  impurity  propa- 
gator. The  dashed  lines  are  the  fermion  propagators  and  the  wavy  and  solid 
lines  are  boson  and  conduction  electrons  respectively. 


the  temperature  dependence  of  the  susceptibilities  agree  with  those  obtained 
from  the  conformal  field  theory  [110,  114]. 

Our  model  includes  a  BCS  pairing  of  electrons  in  the  conduction  band.  In 
general  the  pairing  may  be  either  of  channel  singlet  or  triplet  type.  We  assume 
the  latter  possibility.  The  quantities  we  study  in  this  work  are  not  significantly 
different  for  the  channel  singlet  state.  In  the  limit  of  large  on-site  Coulomb 
repulsion  U  and  for  temperatures  T  <[/ ',  the  model  has  a  simplified  form, 


H  =   E   €AaaCkaa  +  Ef  £  flfa  +  V   £    [4^/4  +  *•* 

k,a,a  a  k,a,a 


+  E  [A(*)cLA,a,-a +  H  +  A  fedf* +E6!fc  -  M 

k,a,o  \   o  a  / 


(5.1) 


The  indices  a  and  a  refer  to  channel  and  spin,  respectively.  The  boson  ba  trans- 
forms according  to  the  conjugate  representation  of  SU(fc).  In  the  limit  A  — >■  oo 
the  unphysical  states  with  the  impurity  occupation  ny  >  1  are  projected  out. 
The  position  of  the  bare  impurity  level  is  assumed  to  be  Et  =  — 0.67D,  and 
T  =  7^(0)1^  =  0.15Z),  where  D  is  half  of  the  band  width.  For  this  parameter 
set  the  Kondo  temperature  is  Tr  =  4.5  x  10_5£)  and  the  impurity  occupation 
number  is  n  t  ~  0.9. 
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The  pairing  correlations  in  the  conduction  band  lead  to  a  nonzero  anoma- 
lous impurity  propagator,  F^a,(r)  =  (TTfa(r)ba(r)fa'(0)bai(0))  (see  Figure 
5.1  ).  This  propagator,  including  the  internal  conduction  electron  line  is  evalu- 
ated self-consistently  at  finite  impurity  concentrations.  Although  Ft  ~  \/N  , 
and  in  NCA  only  the  contributions  0(\/N)  are  retained,  Ft  is  of  order  O(A) 
near  Tc  and  must  be  included  in  the  calculation  of  the  superconducting  transi- 
tion temperature.  There  are  no  anomalous  contributions  to  boson  and  fermion 
self-energies, 

/oo 
def(e)N{e)Gm(u  +  e  +  iO+),  (5.2) 

-oo 


/oo 
de{l  -  f(e))N(e)G0(uj  -  e  +  iO+). 
-00 


jml 


(5.3) 


Here  k  =  2  is  the  number  of  channels.  We  calculate  Ft  in  the  "elastic"  approx- 
imation [56]  in  which  the  internal  anomalous  conduction  electron  propagator, 
see  Figure  5.1  ,  is  evaluated  at  the  external  frequency.  The  slope  of  the  ini- 
tial Tc  suppression  is  obtained  by  including  both  diagonal  and  off-diagonal 
contributions, 


V     dn     7n=0  \     dn     )  x       \     dn     J  2 


n>0 


lmGf{u;n)       Ff(un) 


(5.4) 


The  detailed  form  of  this  formula  in  the  NCA  approximation  for  a  single- 
channel  Kondo  effect  was  derived  in  Ref.      [56].   Here  we  evaluate  analogous 
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expression  for  the  two-channel  model.    The  contribution  associated  with  the 
impurity  spectral  function  is  the  following: 


NQdTc 
dn 
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where  V'  is  the  digamma  function  and  po  and  pm  are  the  slave  boson  and  the 
impurity  fermion  spectral  function,  respectively.  The  second  term  in  equation 
(5.4)  has  the  following  form: 

-  (^fP)     =  fr  jT  <M*/M  -  (m  ^)]/2a(W),  (5.9) 
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(5.10) 
(5.11) 


V(e;  w)  =4  cosh(a;/2r)e-/9eReGo(e)po(e) 

x[e^/2Pm(e  -  u>)ReGm(e  +  w)  -  e~^/2pm(e  +  u;)ReGm(e  -  w)] 

+2sinh(u,/T)e-^[(ReG0(6))2  -  7r2pjj(e)]p»„(e  +  u;)pm(e  -  u). 

(5.12) 

The  expression  (in  <-»  0)  indicates  the  interchange  of  the  fermion  and  boson 
indices,  ?«  and  0,  respectively,  in  the  formula  for  aAu).  The  off-diagonal  con- 
tribution can  be  separated  into  a  spin-flip  and  a  spin-preserving  part,  T^     and 
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Figure  5.2  The  anomalous  part  of  the  conduction  electron  T-matrix  evalu- 
ated at  T  =  4.5xlO-4Tjf.  The  solid  (broken)  line  is  the  non-spin-flip  (spin-flip) 
contribution. 


T2,  respectively.  The  inclusion  of  inelastic  processes  is  not  expected  to  signif- 
icantly affect  the  results  because  the  spin-flip  and  non-spin-flip  contributions 
almost  cancel  [  115,  116]. 

At  low  temperatures  these  off-diagonal  components  of  the  T-matrix  have 
a  peak  at  u>  ~  T  (see  Figure  5.2  ). 

To  evaluate  Tc-suppression  we  use  the  gap  equation, 


=  Vs  /       dujN0{u)—  tanh(w/2Tc), 
Jo  AcJ 


(5.13) 


where  Nq(u>)  is  the  density  of  conduction  electron  states  in  the  normal  metal 
and  Vs  is  the  pairing  potential.  The  full  conduction  electron  Green's  function, 
averaged  over  impurity  positions,  G_1(u;)  =  ojtq  —  6^3  —  A(k)ri,  is  found  from 
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Figure  5.3  The  slope  of  the  initial  Tc-suppression  for  a  superconductor 
with  two-channel  Kondo  impurities.  The  separate  contributions  come  from  the 
diagonal  part  of  the  T-matrix,  7\,  the  off-diagonal  spin-flip,  T^  i  and  the  off- 
diagonal  non-spin-flip,  T2.  The  solid  line  is  the  sum  of  the  three  contributions. 


the  Dyson  equations,  Co  =  u  —  nV^GAuj),  and  A  =  A  +  nV*Ff(u>),  where  n  is 
impurity  concentration. 

As  shown  in  Figure  5.3  ,  the  slope  Aro(0)dTc/dn,  at  n  =  0,  remains  finite 
as  Tk/Tc  —¥  00,  indicating  finite  pairbreaking  in  the  strong  coupling  regime. 
We  include  for  comparison  the  NCA  result  for  the  initial  Tc-suppression  for 
a  superconductor  with  single  channel  Kondo  impurities,  see  Figure  5.4  .  The 
qualitative  form  of  the  Tc  dependence  on  n  in  the  two-channel  case  remains 
practically  unchanged  for  all  values  of  r^-/Tco,  however,  as  illustrated  in  Figure 
5.5  . 

Both  results  are  a  consequence  of  finite  magnetic  moment  of  the  impurity  at 
all  temperatures.  The  maximum  slope  is  reached  for  T^/Tco  — >■  00,  in  contrast 
with  the  single  channel  Kondo  problem  where  maximum  occurs  for  Tj{  ~  Tcq. 
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Figure  5.4  The  slope  of  the  initial  Tc-suppression  for  a  superconductor  with 
single  channel  Kondo  impurities.  The  separate  contributions  come  from  the 
diagonal  part  of  the  T-matrix,  T\,  the  off-diagonal  spin-flip,  T^F ,  and  the  off- 
diagonal  non-spin- flip,  T<i-  The  solid  line  is  the  sum  of  the  three  contributions. 
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Figure  5.5  The  superconducting  transition  temperature  as  a  function  of 
impurity  concentration.  The  inset  shows  the  result  for  a  single  channel  Kondo 
impurity  with  T%  ~  Tc  evaluated  with  Et  =  — 0.67D,  and  T  =  0.15D. 
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Similar  results  for  (dTc/dn)n=Q  were  obtained  in  a  recent  Monte  Carlo  calcula- 
tion [115,116].  If  the  Kondo  scale  is  known  from  other  experiments,  measure- 
ments of  Tc  as  a  function  of  impurity  concentration  may  help  in  identifying 

the  multichannel  behavior.  For  Tcq<T^,  studies  of  the  single  channel  problem 
have  found  a  reentrant  or  exponential  dependence  of  Tc  on  concentration.  If  it 
were  possible  to  vary  the  ratio  T/<-/Tco  in  the  strong  coupling  regime,  increased 
pairbreaking  with  growing  T^/Tcq  would  imply  multichannel  behavior  while 
the  opposite  would  be  true  for  a  single  channel  coupling.  In  unconventional 
superconductors,  for  which  ^2kF(k,uj)  =  0,  both  single  and  multichannel  ex- 
change result  in  qualitatively  similar  form  of  the  Tc-dependence  on  n  in  the 
strong  coupling  regime  and  the  differences  which  may  exist  for  T%  ~  Tcq  may 
be  harder  to  indentify.  Finally,  let  us  note  that  the  breaking  of  the  spin  or 
channel  symmetry  which  takes  the  system  away  from  the  multichannel  fixed 
point  should  reduce  pairbreaking  especially  in  the  strong  coupling  limit.  Given 
the  subtleties  of  this  analysis  we  conclude  that  Tc-suppression  is  not  likely  to 
be  a  sensitive  tool  for  identifying  the  multichannel  exchange. 


CHAPTER  6 
SUMMARY  AND  CONCLUSIONS 

We  have  studied  metals,  mainly  superconductors,  in  the  presence  of  mag- 
netic and  nonmagnetic  impurities.  In  systems  with  power  law  density  of  states 
near  the  Fermi  level  and  antiferromagnetic  exchange  interaction  between  im- 
purity and  conduction  electrons  the  single  impurity  calculation  in  the  large- iV 
mean  field  approximation  [3]  yields  the  normalized  critical  coupling  NqJc  ~  r, 
at  low  r,  and  NqJc  — >■  2,  for  r  — >■  oo.  The  NCA  approach  also  gives  NqJc  ~  r, 
for  r  -C  1.  Recent  NRG  calculations  [33]  confirm  the  result  for  small  r  but 
disagree  at  larger  r.  There  is  an  upper  critical  rc  =  0.5  for  the  orbitally  non- 
degenerate  single  channel  Kondo  problem,  above  which  the  Kondo  effect  does 
not  occur  at  all,  regardless  of  the  strength  of  the  coupling.  The  conclusions 
reached  for  the  same  model  with  nonzero  impurity  concentration  in  the  mean 
field  approximation  are  quite  different:  Jc  =  0  for  any  nonzero  concentration. 

It  would  be  desirable  to  perform  the  NCA  calculations  from  chapter  2  in 
the  Coqblin-Schrieffer  limit,  both  at  finite  temperatures  and  at  T  =  0.  This 
would  allow  a  closer  comparison  of  this  method  with  the  NRG  [33]  results  for 
the  Kondo  model  with  power-law  density  of  states  in  the  conduction  band, 
N(u)  ~  \u\r.  In  addition  to  a  single-impurity  calculation  it  would  be  in- 
teresting to  see  the  self-consistent  NCA  treatment  for  finite  concentration  of 
impurities.  One  would  like  to  know  whether  in  NCA  the  critical  coupling  for 
finite  concentrations  is  Jc  =  0,  as  in  the  mean  field  approach.   Also,  it  would 
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be  of  interest  to  study  the  breakdown  of  the  Fermi  liquid  with  increasing  r  and 
to  use  NCA  to  calculate  dynamical  quantities. 

In  chapter  3  we  have  presented  a  slave  boson  theory  of  Kondo  impurities 
in  superconductors  which  has  the  advantages  of  being  applicable  in  the  Fermi 
liquid  regime  and  being  relatively  easy  to  use  in  calculating  quantities  of  exper- 
imental interest  at  all  temperatures  in  the  superconducting  state.  The  theory 
has  been  shown  to  reproduce  all  of  the  qualitative  features  of  the  physics  found 
by  previous  theories  for  large  T^-,  and  is  asymptotically  in  quantitative  agree- 
ment with  "exact"  NRG  calculations  for  the  s-wave  case.  It  is  furthermore 
capable  of  going  considerably  beyond  currently  available  theories  in  that  prac- 
tical calculations  of  superconducting  response  functions  are  possible  at  low 
temperatures  in  the  superconducting  state,  and  can  be  easily  generalized  to 
the  unconventional  states  of  great  current  interest  in  the  heavy  fermion  and 
high-temperature  superconductivity  problems. 

The  single  exception  to  this  success  is  the  failure  to  properly  describe  the 
s-wave  superconductor  at  very  large  impurity  concentrations,  in  that  the  the- 
ory predicts  an  infinite  critical  concentration.  In  practical  terms  this  is  quite 
academic,  as  all  independent  impurity  analyses  will  break  down  due  to  inter- 
impurity  interaction  effects  long  before  any  putative  critical  concentration  is 
achieved.  Nevertheless,  this  formal  shortcoming  of  the  theory  exists  and  must 
be  addressed.  Preliminary  analysis  of  fluctuations  about  the  saddle  point  has 
convinced  us  that  the  Gaussian  fluctuations  to  the  scattering  amplitude  arising 
from  the  slave  boson  dynamics  will  be  sufficient  to  induce  a  critical  concentra- 
tion, in  analogy  to  the  works  of  Matsuura  et  al.  [52]  and  Sakurai  [53],  and  that 
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the  theory  as  it  stands  is  sufficient  to  describe  the  Fermi  liquid  regime  every- 
where except  for  very  large  concentrations  in  the  s-wave  case.   The  difficulty 
does  not  arise  in  the  unconventional  case. 

In  unconventional  superconductors,  we  have  shown  that  the  phenomeno- 
logical  theories  of  Refs.  [5]  and  [6]  will  be  reproduced  by  the  microscopic 
theory  presented  here  in  the  limit  N  =  2  and  Tj^/Tc  — >  oo.  The  deconden- 
sation of  the  slave  boson  amplitude  prevents  extension  of  the  theory  into  the 
high  temperature  regime,  but  qualitatively  it  is  clear  that  the  effect  of  lowering 
the  Kondo  temperature  is  to  sharpen  the  many-body  resonance  near  the  Fermi 
surface,  but  lower  its  weight.  Thus  the  effect  of  resonant  scattering  leading  to 
low-energy  gapless  effects  in  superconducting  thermodynamic  and  transport 
properties  is  reduced.  Impurities  with  larger  orbital  degeneracy  may  lead  to 
similar  resonances  away  from  the  Fermi  level,  possibly  similar  to  those  observed 
by  Maple  et  al.  [  117]  in  specific  heat  experiments  on  Pr-doped  YBCO. 

The  effect  of  Zn  and  Ni-doping  on  superconducting  properties  of  YBCO 
suggests  that  these  two  dopants  may  be  described  within  a  traditional  Kondo- 
type  picture  in  conjunction  with  the  theory  discussed  here.  Ni,  with  spin 
one,  may  possibly  be  treated  as  a  higher-degeneracy  scatterer  below  its  Kondo 
temperature.  As  we  have  seen,  such  an  impurity  acts  as  a  weak  scatterer 
compared  to  the  N  =  2,  Tr-  3>  Tc  resonant  scattering  case,  a  possible  model 
for  Zn.  On  the  other  hand,  recent  NMR  measurements  appear  to  suggest  that 
a  moment  forms  around  the  Zn  site.  Within  an  isolated  impurity  picture,  this 
would  suggest  that  T%  <C  Tc,  where  the  present  theory  is  not  applicable.  A 
more  plausible  explanation,  however,  is  that  the  local  spin  correlations  induced 
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by  a  missing  Cu  must  be  accounted  for,  as  suggested  by  Poilblanc  et  al.   [98] 
Further  experimental  and  theoretical  work  on  this  problem  is  clearly  essential. 

There  is  by  now  a  considerable  body  of  experimental  data  supporting  the 
picture  of  gapless  superconductivity  in  the  cuprate  high-Tc  materials,  with 
a  residual  density  of  states  and  low-temperature  behavior  varying  qualita- 
tively according  to  the  d-wave  plus  resonant  scattering  model  [16,  118,80,17]. 
This  data  stands  in  apparent  contradiction  to  the  well-known  effect  of  small 
amounts  of  potential  scatterers  on  anisotropic  s-wave  superconductors,  namely 
the  smearing  of  energy  gap  anisotropy.  This  continues  to  hold  even  for  ex- 
tremely anisotropic  systems  with  nodes,  as  illustrated  by  the  simple  theory 
presented  here  for  a  representative  order  parameter.  We  believe  that  this  data 
strongly  suggests  that  the  pairing  is  unconventional  in  these  materials,  but  the 
analysis  presented  in  chapter  4  does  not  as  it  stands  allow  one  to  distinguish 
among  possible  candidate  unconventional  states  (e.g,,  dx2_y2  and  dxz)  without- 
further  quantitative  comparison.  It  should  be  noted  that  time-reversal  break- 
ing unconventional  states  with  a  gap  will  become  gapless  in  the  presence  of 
pure  potential  scattering. 

As  we  have  briefly  discussed  in  chapter  4,  the  major  difficulty  inherent  in 
such  an  analysis  is  the  possibility  that  even  an  apparently  "  inert"  impurity  such 
as  Zn  or  a  vacancy  in  the  Cu-0  planes  may  induce  local  spin  correlations  in 
the  strongly  interacting  electron  system,  leading  to  spin-flip  scattering.  Ruling 
out  gapless  superconductivity  induced  by  magnetic  scattering  then  becomes  a 
quantitative  problem.  Gapless  behavior  in  films  suggests  that  a  resonant  scat- 
tering mechanism  of  some  type  must  be  present  in  order  to  induce  a  significant 
residual  density  of  states  with  comparatively  little  Tc  suppression.    We  have 
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shown,  however,  that  resonant  potential  scattering  does  not  take  place  in  s- 
wave  systems,  and  argued  that  low-energy  resonant  spin  scattering  is  much  less 
likely  than  in  the  isotropic  case.  We  have  furthermore  made  a  crude  estimate 
of  the  importance  of  spin-flip  scattering  in  Zn-doped  YBCO  crystals  which 
indicates  these  materials  are  dominated  by  potential  scattering  and  should 
therefore  exhibit  an  induced  gap  if  the  superconducting  state  is  s-wave. 

We  have  also  shown  that,  although  gapless  behavior  in  thermodynamic 
quantities  qualitatively  similar  to  d-wave  states  is  to  be  expected  for  a  range 
of  impurity  concentrations  in  extended-s  states,  transport  properties  are  quite 
different  in  the  two  states.  In  particular,  residual  T  — >  0  conductivities  (a  and 
k./T)  are  expected  to  be  large  and  to  scale  inversely  with  impurity  concentra- 
tion, in  contrast  to  the  resonant  d-wave  case.  Such  experiments  can  therefore 
be  used  in  conjunction  with  Josephson  measurements  to  settle  the  question  of 
order  parameter  symmetry.  Existing  transport  data  on  YBCO  single  crystals 
appear  to  restrict  possible  s-states  to  those  where  the  average  of  the  order 
parameter  over  the  Fermi  surface  is  nearly  zero.  The  ARPES  experiments 
on  the  BSCCO-2212  material  might  imply  a  strongly  anistropic  s-wave  order 
parameter  with  770  —  0.  The  existing  transport  measurements,  however,  do 
not  support  this  interpretation.  Further  measurements  on  the  latter  system, 
particularly  systematic  doping  studies  to  search  for  an  impurity-induced  gap 
and  to  test  impurity  scaling  of  residual  conductivities  can  help  to  distinguish  s- 
and  d-wave  state.  A  careful  determination  of  the  temperature  dependence  of 
k  over  the  IK  —  30A'  range  in  BSCCO  would  be  of  great  value  in  this  regard. 

In  chapter  5  we  studied  the  effect  of  two-channel  Kondo  impurities  on  su- 
perconducting states.    This  type  of  impurity  is  known  to  have  a  non-Fermi 
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liquid  ground  state  in  a  normal  metal.  Since  the  impurity  never  loses  its  mag- 
netic moment  due  to  overscreening  by  conduction  electrons,  even  in  the  strong 
coupling  regime,  it  is  not  surprising  that  we  find  qualitatively  similar  type  of 
Tc  dependence  on  impurity  concentration  for  all  values  of  Tj^/Tc.  The  ini- 
tial slope  of  Tc  suppression  does  vary  greatly,  however,  reaching  its  maximum 
as  Tx/Tc  — >  oo.  There  is  no  reentrant  behavior  of  Tc  as  a  function  of  n. 
One  should  remember,  however,  that  the  impurity  screening  radius  diverges  as 
T  — >  0  in  the  two-channel  problem.  Therefore,  below  a  certain  temperature 
impurities  cannot  be  treated  as  being  independent  from  one  another.  In  prin- 
ciple the  character  of  the  superconducting  response  in  the  regime  T<Tk  might 
help  to  distinguish  a  two-channel  Kondo  impurity  from  a  single  channel  one, 
provided  we  know  the  symmetry  of  the  order  parameter  of  the  superconducting 
host. 
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